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Abstract 

The elucidation of the properties of the instantons in the topologi¬ 
cally trivial sector has been a long-standing puzzle. Here we claim that 
the properties can be summarized in terms of the geometrical struc¬ 
ture in the configuration space, the valley. The evidence for this claim 
is presented in various ways. The conventional perturbation theory 
and the non-perturbative calculation are unified, and the ambiguity 
of the Borel transform of the perturbation series is removed. A ‘proof’ 
of Bogomolny’s “trick” is presented, which enables us to go beyond 
the dilute-gas approximation. The prediction of the large order be¬ 
havior of the perturbation theory is confirmed by explicit calculations, 
in some cases to the 478-th order. A new type of supersymmetry is 
found as a by-product, and our result is shown to be consistent with 
the non-renormalization theorem. The prediction of the energy levels 
is confirmed with numerical solutions of the Schrodinger equation. 
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1 Introduction 


Elucidation of the properties of the instantons in the topologically trivial 
sector has been a long-standing puzzle. Although there is no doubt as to 
the physical importance of the instantons, it is not trivial to define them in 
this sector. When they are separated by infinite distance they satisfy the 
classical equation. But once the distance between them becomes finite there 
are no classical solutions. Due to the attractive force between the instanton 
and the anti-instanton they can easily collapse to a vacuum. Thus, it is 
impossible to distinguish them from the fluctuations around the vacuum. Of 
course, no pragmatic problem arises in the dilute-gas approximation, but 
the lack of a precise definition has prevented us from going beyond this 
approximation. The main purpose of the present article is to clarify the 
structure of the topologically trivial sector in theories with tunneling, and to 
provide a method which goes beyond the dilute-gas approximation. 

From the early stages of the study of the instanton [[lj, it was evident 
that this difficulty is connected with the fact that the perturbation theory 
in the presence of tunneling is not Borel summable B- ft 0- Both problems 
come from the non-separability of the tunneling and the perturbative effects. 
And if we change the sign of the square of the coupling constant so that the 
force between instantons becomes repulsive, they disappear; the instanton 
configuration is now free from the collapse and at the same time the per¬ 
turbation series becomes Borel summable. Bogomolny || was the first to 
point out that the above “analytic continuation” of the coupling constant is 
the key to going beyond the dilute-gas approximation, and our work gives 
a precise realization of his suggestion in the context of the structures in the 
configuration space. 

In what follows, we will present novel structures for the configuration 
space in the topologically trivial sector. We will perform an explicit anal¬ 
ysis in one-dimensional quantum mechanics, with asymmetric double-well 
potential, but our analysis may apply to quantum field theories with weak 
coupling, though additional complications will be introduced in such cases. 

Our starting point is a geometrical definition of the instantons in the 
topologically trivial sector. As was indicated by Balitsky and Yung ||, the 
instantons in the topologically trivial sector form valleys in the configuration 
space. In Section 2, following Refs. 0-1, 0- ITU] , we review a method to define 
the valleys in the configuration space and present the construction of the val¬ 
ley configurations and their basic constituents, the valley-instantons. Section 
3 is the core of this paper; we will give a detailed analysis of a valley that 
contains a pair consisting of an instanton and an anti-instanton, and show 
that the interplay between the perturbative effects and the tunneling effects 
can be summarized in a unified way in terms of the collective coordinate of 
the valley. The summation over the multi-instantons is carried out in Section 
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Figure 1: The asymmetric double-well potential V(q) defined by Eq.( [2.2| ) for 
e = 5 and g = 0.1. 

4, which leads to the expression for the non-perturbative part of the energy 
spectrum. In Sections 5 and 7, we summarize various tests of our results. 
In Section 5, the prediction of the large order behavior is checked against 
numerical and exact calculation of the perturbative series from the 200-th to 
the 478-th order for a wide range of the parameters of our model. This result 
shows an apparent lack of Borcl singularities at certain values of the param¬ 
eters, which suggests the existence of some non-renormalization theorem. In 
Section 6, we show that the model we analyze has a new type of supersym¬ 
metry, which we dub “A/"-fold supersymmetry”. The non-renormalization 
theorem makes it possible to test the prediction of the energy spectrum. 
These predictions are checked against the numerical values of the energy 
spectrum in Section 7. We conclude with some additional remarks in Section 
8. The four appendices contain some of the supporting calculations. 


2 The valley 


As our object of study we will take a one-dimensional quantum mechanical 
system with an asymmetric double-well potential. We denote the coordinate 
by q, and write the Euclidean action of this system as follows: 


Sfe] 

V(q) 



dr 


\m + vM 


2 \dr / 

9 2 (i - gq) 2 - tgq 


( 2 . 1 ) 

( 2 . 2 ) 


This potential is plotted in Fig.|T] for e = 5 and g = 0.1. For small eg 2 , 
the points q ~ 0 and q ~ 1/g are the local minima of the potential, where 
V(q) ~ 0 and —e, respectively. Thus e is essentially the difference between the 
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energies of the two minima. The parameter g plays the role of the coupling 
constant, since the q 3 and g 4 -terms have coefficients g and g 1 2 , respectively. 
The height of the potential barrier is approximately V (1/2 g) ps 1/(32 g 2 ), and 
throughout this paper we will work in the small coupling regime g ■ Cl. In 
general, the existence of the two minima is limited to the range eg 2 < a/3/18, 
to which we restrict ourselves. We will denote the coordinate of the upper 
and lower local minima by q + and g_, respectively (see Fig.[I]). For the sake 
of definiteness, we assume that g > 0 and e > 0, which is convenient for 
calculation and entails no loss of generality. 

Note that other types of quartic potential for the asymmetric double-well 
type can be made in the form (|2.2|) by a suitable linear transformation on 
the coordinate q and a scale change of the Euclidean time r. Therefore, 
with this form of the potential we are covering a wide range of potentials. 
Furthermore, as we will see later, this form is essential for the argument 
discussing the supersymmetry.^ 

Let us first describe the qualitative features of the valleys. The actual 
construction and calculations will be given in the subsequent sections. 

For e = 0 this model has the instanton and anti-instanton solutions of 
this system, given by the following: 



1 1 

gl + e-C- T d ’ 

(2.3) 


1 1 

g 1 + eC-Tf)' 

(2.4) 


The parameter tj (tj) denotes the positions of the (anti-)instanton. 

In the topologically trivial sector, the simplest valley in the functional 
space of q{r) starts from the vacuum q{r) = 0. On the outskirts of the valley, 
the configuration has an almost-flat region of q ~ 1 /g. When the width R of 
this region becomes infinity, the transient regions become the exact instanton 
and the exact anti-instanton solutions. We will call this valley the II valley. 
Correspondingly, there is also an II valley, which starts from q{r) = 1 /g and 
develops to a pair of / and /. For both valleys the action is a monotonically 
increasing function of R , approaching the sum of the instanton action and 
the anti-instanton action as R —>■ oo. There are more complicated valleys, 
which start either from q{r) = 0 or 1/g and develop to configurations with 
multiple / and /. At the leading order (for small coupling), this part of the 
valley is relevant for the dilute-instanton-gas calculus. 

For e / 0, the valley configurations were found to behave in a qualitatively 
similar manner 10- There are configurations that smoothly connect one 

1 In an earlier paper |J, we have used a (4g 3 g 3 — 3 g 4 q 4 ) term in place of gq in the e 

term of the potential. This was for the sake of the exact relations = 0 and q _ = 1/g. 

But this choice obscures the supersymmetry we present in Section 6 and is thus avoided. 
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minimum to the other minimum, which we call “valley-instanton” and “anti- 
valley- instanton”, and configurations that are similar to the pairs of these 
valley-instantons and anti-valley-instantons exist as solutions of the valley 
equation. For a one-pair configuration with distance R(^> 1) (for preview, 
see Fig.|3]), its action behaves as ~ 2S' < - / ^ — eR— (2 / g 2 )e~ R (Fig.[4|). The 
first term of this action comes from the valley-instanton and the anti-valley- 
instanton, while the second term is the “volume energy” from the region 
between them, where q{r) = g_ (and V{qJ) = —e), and the third term is 
the “interaction” between the instanton and the anti-instanton. When R 
decreases, this configuration reaches the “bounce” solution where the action 
peaks, and then reduces to q{r) = q + for R = 0. At the bounce point, 
the valley line corresponds to the direction of the negative eigenvalue of the 
bounce solution. We will also call this an II valley. There is also an II 
valley (Fig.[5]), in which q{r) = g_ at R = 0, and the action S ( ' II \R ) is a 
monotonically increasing function of R for the whole range of R = 0 ~ oo. 
Specifically, S^ n \R) ~ 2+ eR — (2/g 2 )e~ R for R 1 (Fig.g). 

All these can be seen explicitly and in detail by solving the valley equa¬ 
tion. In the rest of this section, we briefly explain the valley method used, 
and analytically construct the valley configurations. Furthermore, quantities 
needed for later calculations, such as Jacobians and determinants, are also 
obtained. 


2.1 The valley method 


We use the following equation as a precise definition of the valley[[7|: 


dr'- 


v%] is[ 9 ] , is[ 9 ] 


= A- 


5q(r)5q(r') 5q(T') Sq(r) ’ 
This equation can be rewritten in the following form: 

2 


5 


5q(r) 




= 0 . 


(2.5) 


( 2 . 6 ) 


This implies that a solution of the valley equation (|2.5|) extremizes the norm 
of the gradient vector / dT(5S[q\/5q) 2 in the functional subspace of q{r) with 
a fixed value of S, with A playing the role of the Lagrange multiplier. This 
makes this equation suitable for a geometrical definition of the valley (see 

Fig ®- 

More important is the fact that with this definition, the eigenvalue A of 
the second derivative of S is removed from the loop integrations. Initially, 
we separate the integration along the valley line from the whole functional 
integration. After this step, we are left with the integrations over the di¬ 
rections perpendicular to the gradient “vector”, 5S/5q(r), which are the 
loop-integrations. This is done in the following manner. 
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Figure 2: A schematic drawing of a valley in a two-dimensional space. The 
thin curved lines are the contour lines of the action. The thick line denotes 
the valley line. The horizontal arrow denotes the gradient vector G at a point 
q a and the broken vertical line is the integration subspace for tp a , which is 
perpendicular to G. Note that the vector G is not tangential to the valley 
line, which is a general property of the valley equation. 

Firstly, we parameterize the valley by a parameter a and denote the valley 
configuration by g Q (r), and then take a functional integral, 

Z = jJvqe~ s[q \ (2.7) 

where we are adopting the natural unit, Ti = c = 1. (Inserting an operator O 
in the above does not change the following argument.) The normalization of 
the path-integral measure is the same as in Ref. |R| . We define the Faddeev- 
Popov determinant A[</? a ], which is actually the Jacobian for introducing a 
as an integration variable, by the following: 

J dad (^J (p a (r)G(r)dr^ A[<p Q ] = 1, (2.8) 

where <p a {j) = q{r) — q a {r) are the fluctuations over which we will be doing 
loop-integrations, and G(r) is the normalized gradient vector, 


G(r) = / 

By inserting the factor (|2.8|) into the functional integral ( |2.7|) , we obtain the 
following: 

Z = J J da J Vq 5 (^J ip a (r)G(T)dr^ A[ip a \e~ s ^. (2.10) 


5S 




dr 1 . 


(2.9) 
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We expand the action S[q] around </? Q (r) = 0: 


S[q\ = S[q a \ + 


SS[q a ] 

8q a (r) 


ip a (r)dT + 


1 

2 


8 2 S[q a \ 

5q a (r)Sq a (T') 


(Pa(T)(p a (T')dTdT' H-. 


( 2 - 11 ) 

The hrst order term in tp Q {r) vanishes due to the (5-function. The integration 
of the second order term in <p a {r ) yields the determinant of the second func¬ 
tional derivative = 8 2 S[q a \/8q a (r)8q a (T'), in the subspace defined by 

the 8-function. At this order, the Jacobian A[</? a ] is approximated as follows, 


A M= f d j£l G{r)dT = 


dS[q 0 

da 


\ 


dr 


' 8 S[q a 
, Sq(r) 


N A). 


( 2 . 12 ) 


With this, the integral (|2.10|) reduces to, 

da 


Z = J 


V 27t deb D 


Ae~ s[qa] . 


(2.13) 


at the one-loop order. In the above, deb denotes the determinant in the 
subspace described above. Note the apparent reparametrization invariance 
of a, which should exist, as we have not specified the choice of the valley 
parameter a. The factor 1 /a/27t is induced by the fact that the (5-function 
reduces the number of the integrations over q(r) by one. 

The valley equation (|2^ ) dictates that the subspace / < p a {r)G{r)dT = 0 
does not contain the eigenvector of the eigenvalue A. Therefore deb D is 
simply the product of all the eigenvalues of D less A. This means that the 
valley equation ( |2.5| ) converts the eigenvalue A completely to the collective 
coordinate a. This is the key property of the valley equation (|2.5|) , since in 
actual applications one encounters negative, zero, or very small eigenvalues, 
which render the one-loop approximation useless, or at least dangerous. By 
choosing A to be the undesired eigenvalue, we remove it completely from the 
Gaussian (and higher order) integration. 

The valley equation can be made more perspicuous by introducing an 
auxiliary coordinate. First we will rewrite the valley equation (|2.5|) as follows: 


8S Y [q] 

8q(r) 


where S Y is defined by, 




= S[q] 



<*%] 

Sq(r) 


2 


(2.14) 


(2.15) 
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The second term of the above Sy[q] contains a fourth derivative term of q{r) 
with respect to r. It can be removed by introducing an auxiliary variable 
F(t) as follows: 


Sv[?] 



(2.16) 


Taking the functional derivative of Sy[q\ with respect to q{r) and F{r) we 
obtain the following equations: 


dr' 


* 2 S[q] 


dq(r)5q(r 


SS[q] 

5q(r) 

F F (t') = 


F(r), 

A F{t), 


(2.17) 

(2.18) 


which are apparently equivalent to the original valley equation fl2.5p upon 
elimination of F(r). It is evident that any solution of the equation of motion 
is also a solution of the valley equation (|2.18| ) with F(t) = 0. Some of the 
other important properties of the valley method are noted in Appendix A. 

Like any classical solution, the solution of the valley equation (|2.5| ) may 
break some of the symmetries of the system. To restore them, we must intro¬ 
duce additional collective coordinates. The conventional collective coordinate 
method cannot be applied, because, in contrast to the classical solution, the 
existence of the zero modes of the symmetries is not guaranteed for the solu¬ 
tion of the valley equation. We will introduce them in the manner developed 
in Ref. p. In order to achieve definiteness, we will assume that the transla¬ 
tional symmetry is broken by the valley configuration, which is indeed the 
case in the model we treat in this paper. 

Under translation To, the quantum fluctuation </? a (r) = q{r)—q a {r) trans¬ 
forms as 


(Pair) -»• <fa(r) = q(r + r 0 ) - q a (r) 

= <p a (r + r 0 ) + q a (r + r 0 ) - q a (r). 


(2.19) 


This transformation induces an analogue of the gauge transformation on 
S[q a + </? q ]|]T 3[. The basic strategy of our collective coordinate method is 
to introduce the translation r 0 into the path-integral as the “gauge orbit”. 
This procedure automatically guarantees the manifest invariance under the 
translation through the following deformation, and does not require the zero 
mode. 
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To extract the gauge orbit, we will use the Faddeev-Popov technique. 
In place of Eq. ([2.8|) , we introduce the following Faddeev-Popov determinant 
A [ip Q ]: 

J dr 0 j dad (^j ^°(r)G(r)dr^ 6 (^j ^°(r)^(r)dr^ A b«] = 1, (2.20) 

where G{t) is a suitable function independent from r 0 , and gives the gauge 
fixing condition of the translation: 


J Va(j)Q(T)dT = 0. 


( 2 . 21 ) 


By inserting the factor (|2.20|) into the functional integral, we obtain 


Z = J j dr 0 J da J Vqd (^j ^°(r)G(r)dr 
xd (/ ip T ° (t)G {r)dT^j A[(p a ]e~ s[g] . 

As well as the usual gauge symmetry, the following relations hold. 

Aba] = A[^°], S[q a + <f T J] = S[q a + <p a ]. 
Therefore, the integration of r 0 is trivially factored out. 

Z = J J dr 0 J da J Vip a 5 (^j ^ Q (r)G(r) dr 
xd ^I (p a (r)G(T)dT^j A[(p a \e~ s ^ q °‘ +ipa ‘\ 


( 2 . 22 ) 


(2.23) 


(2.24) 


At the one-loop level, we find that the Faddeev-Popov determinant becomes 
as follows: 


AbJ = 


! S? “ W G(r)dr l d jZZg ( r)dr\ 


det 


V 


dr 

dq a (r) 

dr 


da 

G(r)dr ( d jAlG(r)dr 


da 


(2.25) 


Because of the following equation 

dq a (r) f dq a {r) dS 




-dr 


dS 


dr v/ J dr dq a (r) / \ J \dq a (r'), 
it can be simplified as 


dr 1 = 0, (2.26) 


A [ip a \ = 


dq a (t) 

dr 


G(r)dr 


dqjj) 

da 


G{j)dr 


(2.27) 


9 






















Some of the valleys may contain classical solutions, and in the neighbor¬ 
hood of such a solution, the valley has a quasi-zero mode which is different 
from A because of the broken symmetry. To remove this dangerous mode rj 
from the path-integral, we will designate £/(r) as the normalized eigenfunc¬ 
tion of rj: 


dr' „ G(t') = t)Q{t). 


5q(r)5q(r 
With this gauge fixing, we obtain 


Z = J 


27 tV det" D 


A[w|e 


-%a] 


where 


det "D 


det D 
Xrj 


(2.28) 


(2.29) 


(2.30) 


2.2 Valley configurations 


For the action Q2.1|) , the valley equations ([2.171) and (|2.18|) yield the following 
partial differential equations, 


-d 2 q + V'(q) = F 
(-d 2 T + V"(qj)F = A F. 


(2.31) 

(2.32) 


From Eq.( [2.3li) we see that the auxiliary coordinate F acts as an external 
force introduced into the usual equation of motion for q. Eq. (|2.32|) is a self- 
consistent equation for F. 

Since the terms in the action Sy in (|2. 16 ) contain two time-derivatives at 
most, the canonical formalism can be applied. Denoting the “Lagrangian” 
of Sy as Ly, the canonical momentum p q and pp for the coordinates q and 
F are given as follows, 



dLy 

■ 1 rn 

Pi = 

dq 

= q ~X F - 


dLy 

1 

Pf = 

dF 

= -A 9 - 


(2.33) 

(2.34) 


This leads to the following conserved “Hamiltonian” of the system, 


Hv = P q q + PfF - Ly 

= \q 2 - V{q) - i (^F 2 + Fq — FV\q )) . (2.35) 
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The conservation of this quantity can be directly confirmed by the use of the 
valley equations (|2.31| ) and (|2.32|) . 

The valley Hamiltonian Hy in (|2.35|) is useful for numerical calculations 
of the valley configurations. Any II or II valley configuration is symmetric 
around the middle point between the instantons. When solving the differen¬ 
tial equations ([2.31|) and f [2.32|) for a given A, we choose initial values at this 
symmetric point so that the solutions converge as r —> oo. This shooting 
method is somewhat complicated in its operation if the initial data is multi¬ 
dimensional. At the symmetric point, say r = 0, we have g(0) = F(0) = 0. 
The convergence conditions, q{ oo) = q±, F(oo) = 0, are to be satisfied choos¬ 
ing q( 0) and F(0) as appropriate. Therefore at this stage the search has to 
be made in two-dimensional space. The valley Hamiltonian (|2.35| ) reduces 
the degrees of freedom in the initial data by one, since H v = —V(q±) by 
the asymptotic condition for r —> oo. Therefore we only need to choose one 
initial datum, say g(0), to find the valley configuration. Having only one 
parameter to determine, this task is straightforward. 

The numerical calculations were carried out for g — 0.1, e = 5, which 
corresponds to the potential in Fig.[l|. Fig|| is the plot of q a {r) and F a (r) 
of the configurations in the II valley. The values of the action S and the 
eigenvalue A for this valley are plotted in Fig.|], where the horizontal coor¬ 
dinate is 11 q\ | ~ \q- — q+\R for R 1. The configuration drawn by the 
broken line Fig.|| is the bounce solution, which has F(t) = 0. The action 
peaks at this point in the valley as seen in Fig.|j. The corresponding plots 
for the II valley are given in Fig.^]and Fig|| Corresponding to the fact that 
there are no solutions of the equation of motion with the boundary condition 
g(±oo) = g_, the action behaves monotonically, as seen in Fig.^. 

In Fig.|] and Fig|| we readily see that these valleys are qualitatively 
similar to the valleys for the symmetric double-well case e = 0. Particularly, 
for large distances > 1, the valley configurations are in fact made of flat 
inside regions of q = q T and flat outside regions where q = q±. The transient 
regions are of fixed shape, which simply moves outward as R increases (see 
the curves a,b,c of Figs. S0- It is apparent that for R —> oo, each transient 
configuration with boundary conditions q(— oo) = q±,q( oo) = q T becomes 
an independent solution of the valley equation. We call the left transient 
region of the II valley, where q changes from q + to q_, as valley-instanton, 
while the other as the anti-valley-instanton, and vice versa for II valley. The 
behaviour of the eigenvalue A in Fig.[4| and Fig.|j suggests that these valley- 
instantons have an eigenvalue of exactly zero. We note that this property is 
not trivial for e ^ 0, since in that case the valley-instanton is not a solution 
of the equation of motion. In fact, taking the time-derivative of Eq. (|2.31 ), 
we obtain, 

(-% + V"(qj)q = F, (2.36) 
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4(T) 



Figure 3: Configurations in II valley. The symmetric point is chosen to 
be at r = 0, and only the half of the configuration, r > 0, is shown. The 
upper plot shows the behavior of q(r), where the broken line denotes the 
bounce solution. The lower plot is F{r). The labels a, b, and c show the 
correspondence between q{r ) and F(r). 



H 


M 


Figure 4: The action S (upper plot) and the eigenvalue A (lower plot) of 
the II valley. The labels show the action of the corresponding solutions 
in Fig.|3j. The horizontal axis is a valley parameter defined as a = ||g|| = 
fo° \q( t ) ~ q+\dr. For ||g|| 1, ||g|| — \q- — q+\R■ The peak of the action is 

given by the bounce solution shown in Fig.[|. 
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9(T) 



F(x) 


Figure 5: Configurations in II valley. 
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Figure 6: The action S (upper plot) and the eigenvalue A (lower plot) of the 
II valley. There are no solutions of the equation of motion under the bound¬ 
ary condition g(±oo) = g_. Therefore, the action behaves monotonically. 
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which shows that q is not a zero mode. This is as it should be, since simple 
translation of the valley-instanton does change the action by the correspond¬ 
ing volume energy. 

Eq. fl2.36|) , however, can be combined with Eq.( [2.32|) to show that A = 0. 
For this purpose, we multiply q on Eq.( |2.32|) , integrate over r from — oo to oo, 
and carry out partial-integration twice to move d 1 2 to act on q. At r — ±oo, 
q of the valley-instanton is at the extrema q± of the potential. Therefore the 
external force F goes to zero for r —*■ ±oo (see Fig.|3| and Fig.|5|). The surface 
terms of the partial integrations vanish due to this property. We thus obtain 
the following: 

/ FFdr = A / qFdr. (2.37) 


The left-hand side of this equation is evidently zero according to the boundary 
conditions. On the right-hand side, the integral qFdr is not zero, as is 
seen in Fig.|| and in Fig.|5]. Therefore, for the valley-instanton we find that 
A = 0 and the auxiliary coordinate F is the zero mode. Note that this proof 
does not apply to the II and II configurations. This is because for those 
conhgurations the integral qFdr is zero due to the fact that q{r) and 
F(t) are even functions around the symmetric point. And in fact, A is not 
identically zero for the II and II valley. 

Although the global features are evident from the numerical results pre¬ 
sented in this section, we need analytic expressions of the action and other 
related quantities. Of special importance are the interaction terms of the 
actions of the II and II conhgurations. In the following we carry out an 
analytic construction for these conhgurations, using the fact that the valley- 
instanton possesses an exact zero eigenvalue. 


2.3 Valley-instanton 

The existence of the zero eigenvalue A = 0 and the zero mode F enables 
us to carry out the construction of the valley-instanton. We will hrst carry 
out the analytic construction of the valley-instanton. For this purpose it is 
convenient to use the rescaled coordinate q = gq. The action ( |2.1| ) is written 
in terms of this q as follows: 


S[q\ = - 2 -%], S[q] = dr 


1 / dq' 

2 l dr. 


+ V{q) 


(2.38) 


V(q) = V 0 (q) + eg 2 V 1 (q), V 0 (q) = -q 2 (l - q)\ V^q) =-q. (2.39) 


This implies that in the following expansion in e, the actual perturbation 
parameter is eg 2 , which is natural considering that this is a dimensionless 
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quantity. We can also define a rescaled auxiliary coordinate F = gF and 
rewrite the valley equations (|2.31|) and ( |2.32|) as follows: 


~d 2 q + V'(q) = F 
(~d 2 + V"(q)) F = 0, 


(2.40) 

(2.41) 


where we used the fact that A = 0. (In the above the prime denotes the 
derivative with respect to q.) Let us obtain the solution of this set of equa¬ 
tions as perturbative series in eg 2 , 

q = qo + eg 2 qi + {eg 2 ) 2 q 2 + ■■■, (2.42) 

F = eg 2 Fi + (eg 2 ) 2 F 2 + • • •. (2.43) 

As the zeroth order, we have only the following equation: 

- d% + V'(q 0 ) = 0, (2.44) 


whose solution is the (scaled) ordinary instanton q 0 = gqlpir: 0). (We choose 
7/ = 0 for the ease of notation in this analysis. It may be reintroduced by 
replacing r by r — tj at the end.) 

At the first order in eg 2 , Eqs. (l2.40|) and (|2.41|) leads to the following: 


(-a? + L"te„)) = u-Uteo), 

(-0? + h"(«,)) f, = o. 


(2.45) 

(2.46) 


Eq. (|2.46|) dictates that i 7 ) be proportional to the zero mode of q o, which is 
%■ 

F\ = c±q 0 . (2.47) 

The proportionality constant Ci can be fixed by the following consideration: 
For Eq.( p.45| ) to have a solution, the zero mode of the operator — d 2 + V^'lo) 
should not be contained in the right-hand side. Therefore the constant c± 
should be such that Fi cancels the zero mode in V[ (g 0 ). In fact, multiplying q 0 
on both sides of Eq. (|2.45| ), integrating over r from — oo to oo and performing 
partial integration twice, we find that 


/ OO . _ /“OO . _ 

q^dr- / q Q V[{q [} )dr 

-oo J — OO 

= Cl [ q 0 dq 0 - [ Vl(q 0 )dq 0 = \c x + 1, (2.48) 

Jo Jo b 

which Exes ci as C\ = —6. (It is useful to note that % = qo(l — %) in the 
calculation of the first term.) The expression for q\ is obtained by the same 
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procedure as above, but integrating only from — oo to r, which leads to the 
following: 


- f) + Qo = V'(q 0 ) = ho- (2.49) 

The general solution of this equation is hi = l+c 2 <? 0 , where c 2 is an integration 
constant, which remains arbitrary at this order. 

At the second order in eg 2 , we have the following: 

(-dr + K'(ho)) h 2 = h - ^C(ho)h? - W'(ho)hi, (2.50) 

(-^ + K'(ho))F 2 = -C(ho)?iFi. (2.51) 

The parameter c 2 in hi is hxed by the condition that zero mode be absent 
from the right-hand side of Eq. (|2.51|) . Straightforward calculation shows 
that this leads to c 2 = 0. Proceeding in a manner similar to the first-order 
calculation, we find that F 2 satisfies the following equation: 

%F 2 - q 0 F 2 = - 18 ho- ( 2 - 52 ) 

The general solution of this equation is F 2 = (18r + C 3 )h 0 , where the constant 

C 3 is a free parameter. It is straightforward to show that the condition for 

the absence of the zero mode q n on the right-hand side of Eq. (|2.50|) fixes the 

. 2 

parameter C 3 as C 3 = 0 , using the fact that the integral q 0 rdr vanishes 
since the integrand is an odd function of r. 

Summarizing the perturbative solution obtained so far, we have 


q(t) = qo(r) + eg 2 + 0 [(eg 2 ) 2 ) , (2.53) 

F(t) = -6 eg 2 q 0 (r) + 18 (eg 2 ) 2 Tq 0 (r) + O [(eg 2 f) . (2.54) 

The (eg 2 ) 2 term in F indicates that these perturbative series are valid only 
for |t| 1/eg 2 . 

For r —>• ± 00 , we may solve the linearized valley equation, 

[-d 2 + V"(q ± )) Sq = F , (2.55) 

[-d 2 + V"(q±)) F = 0, (2.56) 

where 5q = q — q±, The general solution of these equations are, 

q(r) = 9 , + C^+^reW (2.57) 

F(t) = F T e TU * T , (2.58) 


%<1\ ~ mi = -6 ( y 
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for r —> ±oo, respectively. The coefficients F±, C± are arbitrary constants 


and u± = y V"(q±). 

For small eg 2 , q + = eg 2 + O ((eg 2 ) 2 ) and g_ = 1 + eg 2 + 0((eg 2 ) 2 ). 
Also, u± — It 3eg 2 + O ((eg 2 ) 2 ), which means that the linearized solution 
( |2.57| )- ([2.58p is valid for 1 |r|. Therefore, we find an overlapping region 

1 <C \t\ <C 1 /eg 2 , where both Eqs. Q2.53D , ([2.54]) and Eqs. (|2.57|) , (|2.58|) are 
valid. In this region the exponentials in the asymptotic solutions Eq.( 2.57 ) 
can be expanded in terms of eg 2 r and be matched with Eq. fl2.53|) -~( [2.54|) . This 
matching indeed works out and fixes the constants F±,C± as F± = —6 eg 2 , 
C± = ±1. 

In summary of the analytic construction of the valley-instanton configu¬ 
ration, we have the following: 


q(r) = < 


eg 

1 


c u} + T 


g 

l 


+ 


3 eg 


-re 




CU + 


+ 


glTe~ 

1 1 

-Teg -e 

g g 


eg 


if r -C — 1 


if —1/eg 2 -C r < 1/eg 2 


3eg 

-re 

(U_ 


+ • ■ ■ if r » 1 


(2.59) 


F( t) = 


—6 ege UJ+T -I- 

1 - 3eg 2 r 


-6 eg 


(1 + e -T )(l + e r ) 


-6 ege~ u - T + 


if r -C —1 ; 

+ • • ■ if —1/eg 2 < r < 1/eg 2 ; 

if r » 1 . 


(2.60) 


From this solution we can evaluate the action S. First let us calculate 
the contribution of the solution ( |2.59| ) for |r| -C 1/eg 2 - 

S[q 0 + eg 2 q 1 ] = 5 0 [5o] + eg 2 S 1 [q 0 ] + eg 2 J ^°|^ gi(r)dr -|-, (2.61) 

where we used unsealed quantities qt = qi/g , etc. And is the contribution 
of the Ei-term in Eq. (|2.39|) , while S'o denotes the rest of the terms in the 
action. We find that the second term vanishes. The third term is also zero 
due to Eq. (|2.44|) . In the asymptotic region, |r| ^ 1, the linearized solutions 
should be used, which simply yields the volume energy plus the higher order 
terms in eg 2 . Therefore, the action of the configuration ( |2.59| ) restricted 
in region r G [— T/2,T/2] (T 1) is given by S' ~ (l/6g 2 ) — eT /2 + 
■ • The second term of this expression is the “volume energy” due to the 
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energy difference between q = q + and g_, while the first term is the leading 
contribution to the valley-instanton action S^’^: 


£«.(-0 



1 + o ((t 9 2 ) 2 ) 


(2.62) 


We will now turn to the evaluation of the Jacobian ( |2.12| ). The numerator 
is the squared norm of F(t), the leading contribution to which comes from 
the integration near r ~ 0. Substituting the expression (|2.6(J| ) in this region, 
we find that the leading contribution is, 


/ OO 

F 2 dr 

-OO 


6eV+ •••• 


(2.63) 


Therefore, we find that when we choose the position coordinate of the valley- 
instanton as the valley parameter a, the Jacobian (|2.12|) is given by the 
following: 

£ 1 [l + 0(eg 2 )] . (2,64) 


A = 



F 2 dr 




This leading term is the same as that of the ordinary instanton for e = 0. 

The factor det 7 D can also be evaluated by the use of the fact that A = 0. 
The derivation is given in Appendix B. The result for the valley-instanton 
located at r = tj in the range r G [—T/2,T/2] (T 1) is the following (see 
Eq.(lBJ5|)): 


det'(-d 2 + V ") _ e («--«+)(r/ 2 -rr) 
det(—+ cu+) 2 u_F + F_ 


drF 2 (r). 


(2.65) 


Note that the exponential factor in the ratio of the determinant gives the right 
quantity for the zero-point energy contribution. In fact, the normalization 
factor J of the functional integral ( |2.13|) is given by the following: 


J 


^/det(—<92 +c4) 


_ e -^+ T / 2 'y 


+j 


( 2 . 66 ) 


where Y + is a constant dependent on the definition of Z. For the partition 
function Z = lim^oo Tv(e~ HT ), 


T+ = 1, 

and for the amplitude Z = hiri 7 ^ oc (g + e~ ff7 \q+), 

(2.67) 

T + = |T(0;o; + )| 2 , 

(2.68) 
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where 'I' (q\ u) is the normalized Schrodinger wave function of the ground state 
for the harmonic oscillator of frequency u and therefore \l/(0;u;) = 

Thus we hnd that 


J 


T, 


\Jdet\—d 2 + V") 

where 


uj + (T 


2 V 2 


cn_ (T 


ex P \ —yU + r d-tr o“ T ' > ( 2 - 69 ) 


2 V 2 


«W = 


2cn_Fi F_ J—c 


F 2 dr. 


(2.70) 


For the anti-valley-instanton, det 7 D is evaluated in a similar manner, but 
slightly modified due to the difference of the boundary condition at r = ±oo. 
Eqs.( |2.65|) , (|2.66|) and ( |2.69| ) become as follows, 


det'(-d 2 + V") e (“+-“-)(T/ 2 -T T ) oo 


det(-c>2 +lu 2 _) 2uj + F + F„ 


cItF 2 (t), 


J 


^det(-d 2 + uZ) 


= e -^- T/2 T_, 


J 


T _ 


sjdet\-d 2 + V") ZZ 7 ) 

Here 


exp 


LU- (T 


2 V 2 


+ T r ) - 


u> + (T 


2 V 2 


- r f 


= 


2,U+F+F— J— cx 

and for Z = lim^^oo Tr(e _HT ), Y_ is 

Y_ = 1, 

and for Z = hm T _ s . 00 (g_|e - ' ffT |g_) 

T- = 1^(0; u;_)| 2 . 


F 2 dr, 


(2.71) 

(2.72) 

(2.73) 

(2.74) 

(2.75) 

(2.76) 


The results (|2.69|) and ( |2.73| ) properly represent the zero-point energies 
u±/2 at q = q±. At the leading order in eg 2 , we hnd that u; + = u;_ = 1 and 


— k = 


12 ’ 


(2.77) 


thus the difference between the results ( p.69|) and (|2.73|) disappears. 
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Figure 7: The II valley configuration with a large separation. The symmetric 
point is chosen to be at r = 0, so that 77 = — tj. 


2.4 II valley configurations with large separations 

The rescaled valley equations can be solved for the 11 valley configurations 
in much the same manner as above. The calculation is especially straight¬ 
forward for the leading terms in the eg 2 -expansion, which we will carry out 
in the following.^ 

For this purpose, we will neglect the V\ term in Eq. ( [2.391) . Furthermore, 
we will confine ourselves to the outskirts of the valley, where the eigenvalue 
A <C 1. This allows us to solve the valley equation in a perturbation by A: 


?(t) = Qoir) + Xqi(r) + ■ ■ ■, (2.78) 

F(t) = AFi(r) + X 2 F 2 (t ) + •••. (2.79) 

We will nominate the symmetric point of the configuration as r = 0, and solve 
the equation for the region r > 0 (see Fig.[7]). Therefore we choose % (r) as 
the ordinary anti-instanton solution at r = r/(^> 1); q 0 (r) = gqo\r',Tj). 
The parameter 77 will be related to A through the boundary conditions, 

£(0) =F( 0) = 0, (2.80) 

q(oo) = F( oo) = 0, (2-81) 

for A —> 0 as tj —> cxd. The equations we obtain at each order of A are as 
follows: 


O(A 0 ) : 

-d% + Vo (g 0 ) 

= 0 , 

(2.82) 

0(A X ) : 

(-dl + V''(q 0 j 

) Qi = Fi, 

(2.83) 


(-^t + K ) /, (? o ) / 

)F ± = 0, 

(2.84) 

0(A 2 ) : 

f -< 9 2 + V"(q Q ) 

) F 2 + VfiqofaF! = F v 

(2.85) 


2 This interaction term coincides with the one obtained in Ref. Q. However our deriva¬ 
tion is based on the valley method. 
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The first equation (|2.82| ) is satisfied by the anti-instanton solution. In the fol¬ 
lowing we first solve Eq. (|2.84|) for F \, then Eq. (|2.83|) for q \, and subsequently 
Eq. (|2.85|) for F 2 , satisfying the boundary conditions (|2.81|) at r = oo. At the 


conclusion of the process we require the boundary conditions ( |2.8C| ) at r = 0, 
which read; 


?o( 0) + A ^i(°) = °> 
Fi(0) + AF 2 (0) = 0. 


( 2 . 86 ) 

(2.87) 


The solution of Eq. ([2.84| ) is simply given by the translation mode of the 
anti-instanton q 0 (r)(= — , 0i(r)): 


Fi(t ) = Ci-01 (t). 


( 2 . 88 ) 


This evidently satisfies Fi(oo) = 0. The constant C\ is to be determined by 
the boundary condition at r = 0 later. 

In order to solve Eq. ([2.83|) , it is convenient to define a Green function, 

G(t,t') = i(^i(r)^ 2 (r')0(r-r / ) + - r)), (2.89) 

where '0 2 ('7') is the other zero mode solution of Eq.( 2.84 ), which is chosen so 
that the Wronskian is, 


w = - 01^2 - ^ 1^2 = 2 , 

and therefore has the following asymptotic behaviour: 

( e T ~ Tr for r — tj 1; 


-02 (r) ~ 


-e ( r for t — Tj ■< —1. 


The solution of Eq.( [2.83 ) is then written as follows, 

POO _ _ 

q 1 (r)= dr'G^.T^F^ir'). 

Jo 

It is straightforward to show that gi(oo) = 0. Further, 


Cl 


?i( 0 ) = tt^2(0) / fa( t) dr ~ 


cie J 

~F2 


From this, we find that the boundary condition fl2.86|) leads to, 

12 


ci = —e 

A.9 

Eq. (|2.85|) can be solved as follows: 


-2 TJ 


(2.90) 


(2.91) 


(2.92) 


(2.93) 


(2.94) 


F 2 (t)= dr / G(r,r')5(r / )F 1 (r / ), 

Jo 


(2.95) 
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where B{t) = (1 — 1q (go)gi)- Just as before, this satisfies F 2 (oo) 
t = 0 , we find that 


> 2 (0) = |^2(0 )B, 


0 . At 
(2.96) 


where 


is the quantity of 0 ( 1 ). 
(12.87|) leads to 


roc 

B = / B{T)$l{T)dT, (2.97) 

Jo 

From this, we find that the boundary condition 

2 o 

A = —-e“ 2r h (2.98) 

JD 


The factor B is evaluated to be 1/12 in Appendix C. 

Next we will evaluate the contribution of the region r € [0, 00 ] to the 
action as follows: 


cr 1 cr 1 , % m , ^ , 15s [Qo m , 12 , 

S[q\ = S[q 0 +—-(Agi + •••) + - . 2 (Agi + • • • +••• 

dgo 2 dqo 

A 2 r°° A 

— B[q 0 ] + — / qiFidr — — go(0)gi(0) + • • •. (2.99) 

2 Jo 2 

where we used the fact that the linear term in Agi + • • • vanishes due to the 
equation of motion and the last term is the sum of the surface terms. The 
first term on the left side is given by the following: 


S [go] 




1 

2 


e 


— 2 Tj 


( 2 . 100 ) 


The second term is of order of A 2 ~ e“ 4r/ , and is of next-to-leading order. 
The last term is evaluated as follows, 


A 

2 


9 o( 0 )gi( 0 ) 


A 

2 



^-V ; 2 ( 0 ) f (r) 2 dr ~ 
2 Jo 



( 2 . 101 ) 


Combining all contributions, we find that at the leading order of eg 2 , 

S {IT) = 2 S {1) -eR + S£P, (2.102) 

where is the interaction term, 

slip = (2.103) 

r 

We have denoted the distance between / and I by R = 77 — tj\ = 2 77 . In 
the same way, the II configuration with a large separation R between the 
anti-valley-instanton and the valley-instanton is given by the following at the 
leading order of eg 2 : 

S("> = 2 5 (J) + eR + , (2.104) 
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where S^P is the interaction term, 


(J7) _ R 

°int — 2 C 

g 2 


(2.105) 


Finally, let us evaluate A[yi R ] in Eq. (|2.27|) and det " D in Eq. (|2.30|) . (Here 
we denote (p a as (pn.) We choose the eigenfunction of the second lowest 
eigenvalue of D(t,t') as the gauge fixing function Q{t). This choice enables 
us to be completely free from the difficulties that come from the quasi zero 
mode. The configuration on the II valley constructed above is 


eg 2 ^0 

F (II \t) R =°°6(-t)\ q ( o\T] 7-j) - 6{t)\ q { o\T] rj), 

eg 2 ^0 

where 77 = —Tj and A is given by Eq. (|2.98D . Therefore, we fold 

f dq l {T) G {T)dT «=” * . 

J dR tg 2 ^o -\/ 12 g 2 

The gauge hxing function behaves as follows: 


(2.106) 

(2.107) 


s~t / \ R —^00 

2( r ) 2 = n 

eg 2 ^0 


Consequently, we hnd 




0(-r)qo\r] 17 ) + 0 (r)g^ J (r; rj)) 2 dr 


TO, 


(2.108) 

(2.109) 

( 2 . 110 ) 


dq R 


(il) 


T ) s-> 1 \ 7 R—>oo / 1 

— e ^ dT , s r^3ir 


Combining Eq.( |2.108l ) and Eq.f p.llll ), A [</?#] is evaluated as 

R —>-oo 1 


( 2 . 111 ) 




eg 2 —7-0 6 g 2 


( 2 . 112 ) 


Note that the factorization property at R —> 00 holds: The above A[<^r] 
coincides with the product of the Jacobian ( |2.6Jj ) of the valley-instanton and 
that of the anti-valley-instanton. 

The one-loop determinant can be evaluated from the following equation: 

J 


J 


det (—<92 + u+) 


^det"(—d^ + V"j) yjdet(—d^ + u+) \ ^ et ( + V ii) 


(2.113) 
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where Vjj = V(q^ 1 '). Under the above choice of G(t), the following relation 
holds: 

det"(-d 2 + Vj'j) r-, oo det'(-ff + V[} det'(-ff + Vf) 
det(— d 2 + u\) det(— + uj\) det(— ’ 

where Vj{Vj) is the potential of the (anti-)valley-instanton background. From 
the valley-instanton results (|2.65|) and (p.66|) , we find 


J 


AN fl rjx e 


—T/2 


2n^/det"(-d 2 + Vj'j) ^o 


*9 


(2.115) 


where T = lim^^o Y±. 

For the II valley, we find the following in the similar manner: 


eg 2 -,0 6 g 2 


J 


27 i^jdet"(-d 2 r + Vfj) e 9 


R e~ T ' 2 

A N Tjx^r. 

2 ^o ng z 


(2.116) 

(2.117) 


3 II valley and II valley 


Consider the transition amplitudes between one of the local minima of the 
potential and itself: 


Z 


lim (q+\e HT \q+), 

1 —>oo 

lim (q\e~ HT \q_). 

oo 1 1 ' 


(3.1) 


In general, all configurations made of the same number of / and / are relevant 
for the amplitude, but we will first evaluate this against the background of 
the II and II valleys. Due to the restriction of the background, only a few 
states in the Hilbert space can be treated in an appropriate way, but the 
essence of the physics in the topologically trivial sector already appears in 
this background. The general configurations will be taken into account in 
Section S3 Due to the boundary condition, the II valley is relevant for the 
former amplitude of Eq. (|3.1|) , and the II valley for the latter one. 

As was shown in Section except for the collective coordinates of the 
valleys, the path integral in the background of these valleys can be performed 
with Gaussian integrals. Then we obtain 


Z = lim J 

T— XX) 


= lim IT 

T— > OO 


dr 0 dR 
27 tV det" D 


—S(R) 


X 


a-T/2 

ng 2 


A[ygfj]e 


dR[T - R)e-^ R)/9 \ 


(3.2) 
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where R is the collective coordinate of the valley corresponding to the relative 
distance between the valley-instanton and the anti-valley- instant on ,[] and we 
have abbreviated T(0;u;) with u> = 1 as T(0). (Namely, d'(O) = (1/7r) 1//4 .) 
The factor before the integration of R comes from Eqs. fl2.115Q and (|2.117| ). 
Here we have also performed the integration of the collective coordinate of 
the translational symmetry. Since the configurations in the valleys must be 
one of the minima of the potential at r = ±T/2, a non-trivial factor (T — R) 
results from the integration. As was shown in the previous section, the action 
of the valley behaves as follows: 


S(R) = 
for the II valley, and 
~S(R) = 


(jR 2 at R —> 0 ; 

- 2e~« — eg 2 R at R —> oo, 

3 


aR 2 


at R —> 0 ; 


- — 2e R + eg 2 R at R oo, 
o 


(3.3) 


(3.4) 


for the II valley. Here a is a constant which will be determined later. 

For both of the valleys, the configuration at R = 0 is the vacuum and at 
R ~ oo the well-separated valley-instanton-anti-valley-instanton pair. These 
configurations are smoothly connected along the valleys. In the following, we 
will show that these geometrical structures in the configuration space lead 
to the nontrivial structure of the quantum physics in the topological trivial 
sector. 


3.1 The e = 0 case 


First, let us consider the case of e = 0. We will now carry out the Borel 
transformation of the path-integral, which was first suggested in Ref. ||E) 

El 


and later developed on the basis of the valley method 


Changing the 


integration variable R to t = S(R ), the transition amplitude becomes 


Z= lim |*(0)| 2 X 

x— > 1/3 


e ~ T/2 

Tig 2 


>Cv 


dtF(t)e 


-t/g 2 


(3.5) 


3 Strictly speaking, only when R l$> 1, R coincides with the relative distance between 
the valley-instanton and anti-valley-instanton. For other values of R, we define R so as to 


satisfy Eqs.(3.2)-(3.4). 
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where x = S(T ) = 1/3 — 2e T , Cy = [0, x]. From Eq. ( p.3| ) , it is found that 
F(t) behaves as follows: 


m 


' T 
a/4 at 

1/3 — t \l/3 -x) 


for t —> 0 ; 
for t —> 1/3. 


(3.6) 


To obtain the full form of F(t), a detailed analysis of the valley is needed. 
Instead, we will simply assume that the form of F(t) is as follows: 


m 


m 

a/4ct£(1/3 — t ) 


In 


/ 1/3 - t \ 
\l/3-x) 


where f(t ) is a function that satisfies 


/( 0 ) 


T 

3 ln(l — 3x) ’ 



(3.7) 


(3.8) 


The integral (|3.5| ) contains both the perturbative contribution at t ~ 0 
and the non-perturbative one at t ~ 1/3. To separate the perturbative and 
non-perturbative contributions, we can deform the contour Cy to the sum of 
Cp and Cnp as is shown in Fig.||. Then the amplitude becomes 

—T/2 

+ lim ^(O)! 2 x -—— / dtF(t)e- t/9 \ (3.9) 

x->i/3 ng z JC N p 


Note that there is a significant resemblance between the first term on the 
right-hand side and the formal Borel summation of the perturbation se¬ 
ries. Therefore, we identify the first term of Eq. ( | 3. 9j ) as the formal Borel- 
summation of the perturbative series, and the second term as the non- 
perturbative contribution. As we shall see, this decomposition naturally ex¬ 
plains the interplay between them. For simplicity, we denote the Erst term 
as Zp and the second as Z^p- 

According to the above assumption, the Erst term of Eq. (|3.9|) is given by 
the ordinary perturbation theory. Schematically, it is given by 

Z P = lim |T(0)| 2 x e- T / 2 (l + P 0 (g 2 ) + P 1 (g 2 )T + P 2 (g 2 )T 2 + • • •) ,(3.10) 

1 —KX> v / 


where Pi(g 2 ) (i = 0,1, 2, • • ■) are polynomials starting with g 2 -terms. In order 
to avoid contradicting the assumption, the parameter a is determined as 


1 

4t ig 2 T 2 ' 


(3.11) 
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Figure 8 : Deformation of the contour Cy to the sum of Cp and Cnp. 
The second term of Eq.( |3.9| ) can be evaluated as follows, 

—T/2 rx 


Z N p = lim |\E f ( 0)|' 2 x 
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(3.12) 


where s = t — 1/3, and y is the angle of the path Cnp at infinity. In the 
above we have retained only the leading contribution in g 2 . Therefore, we 
obtain the following equation. 


Z = Z, 


J NP 


= lim 

rp _. „ Z_ J 


l *( 0)| 3 


2 


i + Po(g 2 ) 


+aM^+7ln(-l)+5{ln(-l 
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1 + T < Pi(g 2 ) + — a 2 < 7 + In ( —— ] + • 


+^-{Ato 2 ) + « 2 {i + •••}} + 


(3.13) 


where a = e' 1 ^ 2 /yjng 2 and Pi(g 2 ) = Pi(g 2 )/( 1 + Po^ 2 )) (* = 1,2). 

To obtain the non-perturbative corrections, we can compare this with the 
following: 


z = r lim (?+|e 19+) 

1 —KX) 

= T li m E |^(9+)| 2 e-^ {e=0 ’ JV+=0)T + ---, 


(3.14) 


C=±l 


where F f (e = 0, N + = 0),? = ±1 are the energies of the first two lowest 
states and v h f (g + ) are the values of the normalized wave functions of the 
states at q — q + . These states are perturbatively degenerated as, 


£ t (e = 0,lV + = 0) = - + 0( 9 2 ), 


(3.15) 


and the degeneracy is resolved by the tunneling effect. The non-perturbative 
corrections for these energies are thus found as 


7,NP 


[e = 0, N + = 0) = qa + a 2 \ 7 + In —- 


r 


(3.16) 


where = ±1. In addition, we obtain the non-perturbative correction for 
the values of the wave functions at q = q + : 


I^c(9+)Inp 

a 2 , l2 I T"(l) 
= — 1 ^( 0)| 2 ’ 


+7ln r ? + H ln (- 7 , 


(3.17) 


Note that there are imaginary terms in Eqs. (|3.16|) and (|3.17|) . As will be 
explained in Section |3.3| , this is because that the perturbative corrections for 
them are non-Borel-summable. 


3.2 The e 7 ^ 0 case 
3.2.1 II valley 

As was seen in Fig.|4|, the action of the II valley S(R) is not an invertible 
function of R. Then, if we change the integral variable R to t = S(R) in 


Eq. (|3.2|) , there arises an apparent singularity. To avoid this, we add eg 2 Si(R) 

eg‘Si = fjj‘ [ d,Tq R (T) 


(3.18) 












from the action and define the new integral variable t as 

t = S(R ) + eg 2 Si(R) = S 0 (R). (3.19) 

Here c[r{t) is the II valley with R. The function Sq(R) is invertible, so this 
causes no apparent singularity. When e —> 0, this definition of t smoothly 
coincides with that in the case of e = 0, since eg 2 S\{R) is zero when e = 0. 

The remaining analysis is an extension of the case of e = 0. Using t, the 
transition amplitude is rewritten as 

Z= lim |T(0)| 2 x -—— f dtF{t)e~ t/9 \ (3.20) 

x >1/3 7T^ J C\r 


where x = Sq(T) = 1/3 — 2e T and Cy = [0, x]. As Sq(R) behaves as 


aR 2 


for R —> 0 ; 


So(R) = { i 


2e R for R —> oo, 


then F(t) satisfies 


T 


Fit) = 


y'Acrt 

1 /1/3 — t N 

l 1/3 — t V 2 / 


for t —> 0 ; 

In ( for t —> 1/3. 


(3.21) 


(3.22) 


.1/3 — x 


We assume the form of F(t) to be 


F(t) = 


fit) 


1/3 — t' 


\/4(jt(l/3 — t)\ 2 

where f(t) is a function which satisfies 

6 ~ e T 


In 


1/3 -t' 
1/3-x 


m = 


3 ln(l — 3x) ’ 


/(i/3) = RFi. 


(3.23) 

(3.24) 


We can identify the perturbative contribution and non-perturbative contribu¬ 
tion in the same manner of the case of e = 0 , and we find the non-perturbative 
contribution to be as follows. 
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(3.25) 
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(3.26) 


where ijj(z) = dlogT(z)/dz is the psi function. 


The result (|3.25|) leads to the following for the non-perturbative contribution 
to the perturbative ground state at q = q + , 


4 + P ) (e,iV + = 0) = 


-a 


T(-e). 


(3.27) 


Furthermore, the non-perturbative contribution to the normalized wave func¬ 
tion of this state is found to be the following 

I'I'(9+)Inp = -a 2 |<I>(0)| 2 

x|ln('~V(- e )-r'(- £ )|. (3.28) 

The case where e=A/">0,A/*GZ requires further discussion, which will be 
given in the next section. 


3.2.2 II valley 

The evaluation of the amplitude in the background of the II valley can be 
performed in the same manner as that of the II valley. The non-perturbative 
part of the amplitude is the following: 


J NP 


(e-l/2)T—1/3 9 2 

lim |T(0)| 2 x---- 

T—>oo 7ig z 



r (e)T 




r(e) 
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(3.29) 
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+ E 

N+= 0 


( 2 T 

e -(N++e)T 

U 2 J 

N + \(N + + e) 2 


From the above equation, we find the following non-perturbative contribu¬ 
tions for the perturbative ground state at q = g_: 




(3.30) 


I*(?-)Inp = -o 2 |*(0)| 2 (-^) 

x jlid-^W)~r'(e)J. (3.31) 

3.3 Singularity of the Borel plane 

An immediate consequence of our decomposition of the perturbative and 
non-perturbative contribution is 


ImZp + IiiiZnp = 0. 


(3.32) 


This is because Z = Zp+Z-^p is real. Thus, this simple equation explains why 
the imaginary part of the formal Borel summation of the perturbation series 
is canceled by that of the non-perturbative contribution. At the same time, 
it also shows that the nonzero imaginary part of the non-perturbative contri¬ 
bution is a necessary and sufficient condition for the non-Borel-summability 
of the perturbative expansion. 

Furthermore, assuming that fit) in Eqs.(JT7) and ( |3.23| ) has the following 
form: 


fit) = f(t,g 2 /t), (3.33) 

where /(t, s ) is an analytic function of t for small s, we can predict the large 
order behavior of the perturbative contribution. To this end, we examine 
the analyticity of Zp{g 2 ) in the complex g 2 -plane. Let 9 denote the phase of 
g 2 , g 2 = \g 2 \e ie . When the phase of g 2 changes to 27T, the contour of Zp(g 2 ) 
changes as Fig.|9. Thus we obtain 

Zp{g 2 e 2m ) = Zp(g 2 ) + lim |T(0)| 2 x e ~— [ (3.34) 

z^i/3 ng z Jc 

where the contour C is given in Fig.[T0|. Therefore, Zp{g 2 ) has a cut on the 
real axis in the complex g 2 -plane. This is the only singularity near the origin. 
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Figure 9: The change of the contour Cp(6) of Zp{\g 2 \e ie ) as 9 is changed from 
zero to 2tc. 


Thus the dispersion relation becomes 
Zp{g 2 


= J_x dz z ^) 


2ni Jc g 2 z — g 2 

1 r°° , ImZp(z) 

- / dz -+ ... 

TT Jo z — g z 

1 f°° , ImZ NP (z) 
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TC Jo 

oo 


z-g* 


+ 


-- 


K n=0 J ° 


(3.35) 


71 =U 

where C g 2 is the contour that circles around z = g 2 , and we have neglected 
the contribution from the singularity far from the origin. Thus, the following 
large order behavior of Zp(g 2 ) is predicted: 


OO 


m 


Zp (g 2 ) = c m g 2 

m =1 

m—*oo 1 / 30 , 2 In^-^NP (9 ) 

c m = -/ dg 


7T Jo 


9 


,2m+2 


(3.36) 


Rewriting this in terms of the energies and the wave functions, we find the 
perturbative contribution to the energy Ep and to the wave function |\E r (<j l ±)|p 
as follows, 


OO 

Ep (g 2 ) = Y, e ™9 
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m—>00 


1 , 2 ImT NP (g 2 ) 
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dg 
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(3.37) 
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Figure 10: The contour of integral 


and 


l^(?±)|p = 

m =0 


Wr 


dg 


iIm|^(g±)|NP 


^ 2771+2 


7T j o g z 

Detailed tests of these predictions are performed in Section |]. 


(3.38) 


3.4 Bogomolny’s trick 

In the early eighties, Bogomolny || suggested a procedure to calculate the 
tunneling effect beyond the dilute-gas approximation. He suggested that by 
the formal analytic continuation of the coupling constant which changes the 
force between instantons from attractive to repulsive, the instantons in the 
topologically trivial sector become meaningful and their contributions can be 
separated from the perturbation theory. The same trick was used by Zinn- 
Justin [|T4[ to obtain the energy levels of the excited states. Our identification 
of the non-perturbative contribution is a precise realization of Bogomolny’s 
suggestion, and it provides an explicit justification for this trick. 

To see this, let us examine Znp((? 2 ) in the complex g 2 -plane. As is shown 
in Fig.|TT|, if we perform the analytic continuation of Z^p{\g 2 \e ld ) from 9 = 0 
to 9 = 7T, the contour for Z N p(\g 2 \e ie ) changes from C'np(O) to C'np(vt). And 
in the weak coupling limit the integral of Cnp(vt) can be well-approximated 
by that of Cy, because the dominant contribution of the integral comes 
from t ~ x. Therefore, in the case of g 2 = \g 2 \e m when the interaction 
between instantons is repulsive, Z^p{\g 2 \e l ' K ) coincides approximately with 
what Bogomolny suggested as a method of evaluation: 


ZMWVn = Z(\g 2 \en 


(3.39) 
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Figure 11: The change of the contour Cnp(0) of Z^p(\g 2 \e ie ) as 9 is changed 
from zero to n. 

Note that the coincidence Z NP (g 2 ) = Z(g 2 ) does not hold for 6 = 0; while 
Z N p(g 2 ) could be complex, Z(g 2 ) never possesses the imaginary part in this 
case. This is because a so-called Stokes phenomena occurs for Z(g 2 ); the 
terms in Z{g 2 ) which are negligible for 6 = 7r become relevant for 9 = 0. Since 
the imaginary part of Z^p(g 2 ) is necessary to cancel the ambiguity of the 
Borel transform of the perturbation theory, the formal analytic continuation 
of Z(g 2 ) which Bogomolny originally suggested should be replaced with the 
real analytic continuation of Z^p(g 2 ). 

One can easily find that Z^p(g 2 ) has the same asymptotic expansion for 
9 = 0 and 9 = n. Therefore, the asymptotic expansion of Z^p{g 2 ) with 9 = 0 
can be obtained by the analytic continuation of the asymptotic expansion of 
Z(g 2 ) with 9 = n. We will use this in the next section. 


4 The multi valley 


4.1 Multi valley calculus 

In this section we will evaluate the partition function Z = linip^oo Tr(e -HT ), 
by summing over those configurations made of several valley-instantons, by 
utilizing the knowledge of the valley-instantons and the interactions among 
them in Section 2, and applying the analytic continuation discussed in Section 
371. This enables us to evaluate non-perturbative contributions to excited 


states as well as the ground state. 

We will take a valley made of n-pairs of the valley-instanton and the 
anti-valley-instanton with periodic boundary condition. Since we perform 
the calculation for g 2 = \g 2 \e in < 0, the force between the valley-instanton 
and the anti-valley-instanton becomes repulsive. Therefore, the configura- 
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Figure 12: The collective coordinates Ri and Ri for n-pair configuration of 
valley-instanton and anti-valley-instanton. 


tions with large separations between valley-instantons dominate. From the 
calculation of the action for one pair in Section 2, we find that the action of 
this n-pair configuration is given by the following: 


S 


n 

3? 


i =1 


2 

9 2 


E e ' Rl 

i=l 




(4.1) 


where Ri is the distance between the i-th valley-instanton and the i-th anti- 
valley-instanton and Ri is the distance between the i-th anti-valley-instanton 
and the (i + l)-th valley-instanton (mod n ). (See Fig.|l^.) The expression 
(|4. 1|) is valid when all the valley-instanton and the anti-valley-instantons are 
well-separated, that is, when Ri, Ri 1. 

We will write the sum of the contributions of the n-pairs of valley- 
instantons to the partition function Z for n = 1 ~ oo as follows, 


Z NP = lim £ « 2n J n , 

T—>oo ' 


(4.2) 


n=1 


where the constant a contains the contributions of the Jacobian ( 2.64 ) and 
the ^-independent part of the determinant ( |2.77| ), and is at the leading order 
in eg 2 , 
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The term ./„ is dehned as follows: 
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(4.4) 
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Eq. ( [4.2|) can be evaluated in a manner similar to that performed by Zinn- 
Justin for the ordinary instanton case [[14[]. First, we rewrite the 5-function 
as follows: 


5 XXiii + £0 - T = 


f'lOO—T] 


\i=l 


2lli J—ioo—r 


dsexp l-sT + s^2(Ri + Ri)\ , (4.5) 


i =1 


where 7] is a positive real number. This allows factorization of the integrals 
over R{ and R t in the following manner. 

| rioo—T] 


J m. 


27TX J — 200 —Tj 


ds K + (s) n K_(s ) n , 


(4.6) 


where K±{s) are, 


K±(s) = / dR exp s±R-\—-e 


-R 


s± 
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S_|_ = s 


r(-»±), 


s- = s + e 


(4.7) 

(4.8) 


This leads to the following expression. 
Znp = 

y rioo—T} 


— r” dse~ Ts y 

2ni J—ioo—iq n 


2ni J—ioo—r] 

By partial integration we find that 


/ ZOO—7] 
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where 


4>(s) = 1 — a 2 K + (s)K_(s). 


(4.9) 


(4.10) 


(4.11) 


From Eq.(4.7), we see that 0(s) has both poles and zeros in s. Denoting the 
poles by s = E 0 (N ) and the zeros by E(N), we find that 


^np — E 


d -E(N)T 


N 


E 

N 


,-E 0 {N)T 


(4.12) 


The poles of </>(s) are given by the poles of the T-functions, s = -Eq (e, N±) 
dehned as follows, 


E { 0 + \e,N + ) = - + JV+, 

E b, JV_) — —e+- + N_, 


(4.13) 

(4.14) 
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which are the perturbative zeroth order energies of the iV±-th excited states 
around q = q±, respectively. The subtraction in Eq. Q4.12D corresponds to 
the fact that our Znp does not contain this zeroth order term (and in fact 
vanishes for g — 0). Therefore the full partition function has only the first 
term of Eq. (|4.12| ). This means that the non-perturbative contributions to 
the energy spectrum are given by the zeros of 4>(s): 


4i(s) = 0, 


(4.15) 


This equation is identical to the result obtained by the WKB approximation 
(this derivation is given in Appendix D). 

Let us solve Eq. (|4.15|) in a perturbation of a -C 1. In the case when 
the poles of the T-functions of K + and K_ do not coincide, the solutions 
s = E^^e, N ±) + jp (e, N±) are found by a series expansion in a 2 : 


E$(e,N ± )=a 


,(— 1)" ±+1 


\ ±e+2N± 


N. 




9 


r(=Fe- N±) +0(« 4 ). (4.16) 


The result Q4.16|) may diverge depending on e and N±. In fact, the poles 
occur when the argument of the T-functions is zero or an negative integer. 
This happens when the zeroth order energies of the states around q + and 
g_ are degenerate, which is allowed when e = A f (A f G Z). In fact, these 
divergences in i%p (e, N ±) are caused by the confluence of the corresponding 
poles of the T functions in K + and Jv_. In cases where this happens, we 
need to use a different expansion scheme in a, E(N) = E 0 (N) + api + a 2 p 2 + 
0(a 3 ). The zeroth order term E(N) is chosen so that the arguments of the T 
functions in K± are zero or a negative integer —N± at the same time, which 
means 

(4.17) 


E„(N ± ) = - + N + = —e + - + AC. 


This situation is illustrated in Fig.|I^. Straightforward calculation yields the 
following non-perturbative correction E ? N p(e, A^ + ) to the energy levels: 


-Ey N p(e, Af + ) = qa< 


\ N + \N_\ 


( 2 _\ n++n - 

U 2 / 


a 2 1 
T N + \NJ. 


f2\ N++N ~ 

U 2 / 



+ 1) — i/)(N- + 1) , 
(4.18) 


where q = ±1 and N_ = N + + e = Ah|_ + A f. The plus and the minus signs in 
the expression Q4.18p correspond to two (approximate) linear combinations of 
the perturbative states in the left and right well with the same zeroth order 
energy. This situation is analogous to the lifting of the degeneracy by the 
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Figure 13: The potential V(q) and the zeroth order eigenstates when 
Eq. (|4.17|) is satished. 


instanton contribution for the symmetric double-well potential. (In fact, in 
the limit e —> 0 the ordinary instanton result for the energy splitting between 
the ground state and the first excited state is recovered.) 

The non-perturbative contributions to the energy eigenvalue obtained 
here contain both a real part and an imaginary part. As discussed in the 
previous section the imaginary part is related to the large order behaviour 
of the perturbative series. We will confirm this in the following section by 
directly calculating the perturbative coefficients and comparing them with 
the prediction of the valley method. The real part of the energy needs fur¬ 
ther consideration: Normally the non-perturbative contribution is obscured 
by the perturbative contribution. The quantities which receive only a non- 
perturbative contribution provide a meaningful test of the non-perturbative 
effect, while other quantities do not provide such a test. We will discuss such 
a non-renormalization theorem in Section 6, and then using that result we 
will test the real part of the energy expression in Section 7. 


4.2 II valley and II valley revisited 


Before proceeding to the various tests of the result in Section [47T| , it is instruc¬ 
tive to go back to the // and // valleys and compare the result of Section 
|3] with that of Section |Q| . If we compare Eq.( p.l6| ) with Eq.( |4.18|) , and 
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with Eq.(|4.16|), it will be seen that there is a complete 


Eqs.( p.27| ) and 
agreement between them. 

In Section [| we did not evaluate the non-perturbative contribution to 
the perturbative ground state at q = q + in the case of e = Af > 0, Af G 
Z. The reason is that in this case the ground state is degenerate with the 
perturbative Af-th excited state at q = q_, thus the mixing of the states 
should be taken into account. Without knowledge of the mixing, the non- 
perturbative correction to the energy levels cannot be uniquely determined.^ 


We, however, find that the result (|3.26|) is consistent with the result 


(|4.18|) if the non-perturbative contribution of the perturbative ground state 
at q = q+ satisfies the following relations: 


-^?,Np(e, N + — 0) — 

. A r r 

TM l ~g* 



I^(0+)Inp =*« 


~\~Ot 


V’(l) - i’i-hf + 1) , (4.19) 


(^(1) - i>{Af + 1)) 


9 \ A/* A/* I 

?)' w+1 >-£5*£h 


(4.20) 


where <; = ±1. We also note that the above non-perturbative corrections 
for the energy and the wave function are reproduced by the following 2x2 
effective Hamiltonian: 


Hnp — 



where 


W = a 



(4.21) 

(4.22) 

(4.23) 

(4.24) 


4 In the case of A/"=0, Z 2 symmetry exists in the system, which makes it possible to 


determine the mixing. This is why we can obtain Eqs.(3.16) and (3.17). 
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The (1,1) component W is interpreted as the matrix element between the 
perturbative ground state at q = q + and itself, and the (2,2) component X 
as that between the perturbative Af-th excited state at q — q _and itself. The 
off-diagonal component Y is the matrix element between the ground state 
and the Af-th excited state. This result supports Eqs. Ql.l^ ) and (|4.20| ). 

Finally, we would like to comment on the non-perturbative corrections 
for the normalized wave function of the ground state, Eqs.( |3.17|) , (|3.28|) , 
(|3.31|) and (|4.20|) . These corrections cannot be evaluated by the calculation 
in Section |4. 1| , where we have considered Z = linnr^oo Tv(e~ HT ), since no 
information concerning the wave function is contained in this partition func¬ 
tion. (See Eq.( [f.l2| ).) According to the argument of Section |T3], we can 
predict that the perturbation series in g 2 for the normalized wave function 
of the ground state should be non-Borel-summable. 


5 Large order behavior of the perturbative 
series 


The expressions of the non-perturbative contribution to the energy levels 
(|4.16| ) and (|4.18|) contain imaginary parts starting at the order a 2 , and are 
given by the following for the non-degenerate case: 


Im 




2 \ ±e+2N 7r 

'g 2 ) T(1 ±e + V)’ 


(5.1) 


where we have omitted the subscripts ± of N for brevity. For the degenerate 
case e = Af, we have 


Im(£ ? , N p(AT,iV)) = -a 2 


7T 


A) 


2N+M 


N\(N+Af)\ 

We note that the result (|5.1|) reduces to ( |5.2|) for e = Af: 

Ill! (EtpW, N )) = Im (£fe>(JV, N + Mj) = Im (B t , NP (V, N)) 


(5.2) 

(5.3) 


The imaginary part of the energy is continuous in e at e = Af. 

We recall the relation between the imaginary part of the non-perturbative 
contribution and the large order behaviour of the perturbative series ex¬ 
plained in Section 3.4. The following relation holds: 


E^\e,N,m) 


1 
7T 



d,j‘ 


M-eSAaQ) 

g2m+2 


(5.4) 
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for the perturbative coefficients of the energy levels defined by 


4 ±} (e,iV)= Y i ^ ±) (e,N,m)g 2m . 


<±), 


(5.5) 


m =0 


Substituting the result ( [5 . 1|) into Eq.( fT4l) , we obtain the following for the 
perturbative coefficient N,m), 


E^\e,N,m) = A (±) (e, N)3 m T(± e + 2N + m + 1) 

1 


x 


1 + 0 


m 


(5.6) 


where the coefficient A^{e, N) is dehned as follows, 


A (± \e,N ) = 


q±e+2N 


it N\T(±e + l + N) 


(5.7) 


We note that the higher order corrections of 0(eg 2 ) and 0(g 2 ) that we have 
ignored in Eq. (|5.1|) contribute only to the 0(l/m) correction, or more specif¬ 
ically the 0(l/(±e + 2N + m )) correction in the result ([TT])f| 

The expression for Ep \e, 1) coincides with the expression obtained in 
Ref. |l^, 19]. In order to conhrm our prediction, we have independently car¬ 


ried out the numerical and exact calculation of the perturbative coefficients 
4 ± ' ) (+ N, m) for much wider parameter ranges and to much higher orders, 
using the methods described in Ref. |^0], [l|, [HJ. Specifically, the method 
given in Ref.JT 8 | works for any excited state by replacing the zeroth order 
wave function by the n-th order harmonic oscillator wave function, and is 
fastest in exact computer calculations. We have calculated the following five 
categories of the perturbative coefficients to the 200 -th order: 


(a) N — 0 (—) level (the perturbative ground state at q = g_) for e = 0 to 
10 with Ae = 0.2 interval. Floating point calculation. 

(b) N — 0 (+) level (the perturbative ground state at q = q + ) for e = 0 to 
20 with Ae = 0.2. Floating point calculation. 

(c) N = 1,2,3 (—) levels for e = 0 to 10.5 with Ae = 0.5. Exact calcula¬ 
tion, 

(d) N = 1,2,3 (+) levels for e = 0 to 20 with Ae = 1.5. Exact calculation, 

(e) N = 4, 5, 6 (—) levels for e = 2.5. Floating point calculation. 

5 An attempt to calculate this type of correction is given in Ref. JlTj. 
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Figure 14: The plot of the ratio of successive perturbative coefficients r m 
(dots) and the fitted curve (solid line) for the ground state at e = 9.8. The 
fitting was done for m = 150 ~ 200. 


In the above, “exact” calculation indicates an algebraic calculation carried 
out in Mathematica. 

We can compare these results with Eqs.( |5.6|) and ( |5.7|) with the follow¬ 
ing procedures. First, in order to check the leading m-dependent terms in 
Eq. (|5.6|) , we take their ratio. 


Ep ±] (e, N, m) 
Ep Z \e, N,m — 1) 


3(± e + 2 N + m ). 


(5.8) 


We use the result of the corresponding perturbative coefficients E m for m = 
150 to 200, and fit them as follows: 


E„ 


r m = 


(3-2 a 3 


a 4 


Em —1 


— clq + dim H-1- — H- 


m 


m~ 




(5.9) 


for all the categories (a)-(e). A plot of r m for the ground state for e = 9.8 
is given in Fig.[L4| as a sample. The fit works to good accuracy for values 
from m = 150 to 200 as can be seen in this figure. The same is true for all 
the cases in the above five categories (a)-(e) we have examined so far. The 
resulting coefficients are compared with their theoretical predictions: 


a 0 = 3(±e + 2N), a 4 = 3. 


(5.10) 


In all the cases, the theoretical prediction ( |5.10|) is conhrmed with good 
accuracy, the maximum error being 2 %. 
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Next, we examine A^\e,N) defined by Eq.( |5.6| ). For this purpose we 
use the following quantity b m : 


b m — (ie + 2N + ?7i + 2)d m _|_i — (ie N m 1)^ 


where 


dm. — 


E„ 


(5.11) 


(5.12) 


3"T(± e + 2N + m) 

This quantity is defined so that it is free from the correction term of order 
l/(e + 2 N + m + 1) to E m . (This is not particularly crucial, since the order 
of our calculation is high.) Thus the numerical fit can be performed with the 
following formula: 


b m — c o + 


c 2 


+ ••• + 


C 8 


m* m 

and is compared with the theoretical prediction: 


8 ' 


Co 


= yl (±) (t,JV). 


(5.13) 


(5.14) 


The result for case (a) is plotted in Fig. 15. The difference between Eq.(5.7) 
and the calculated value is at most 0.1 %, which happens at e = 9.8. For 
case (b), plotted in Fig.[L6], the error is at most 15 % (at e = 20). Cases (c) 
and (d) are plotted in Fig.[L7]. The maximum error is 65 % for the case of 
(20.5, 3). Since this error is somewhat large, we have carried out a yet 
higher order calculation for this case. The error decreases to 32 % at the 
364-th order, and to 17 % at the 478-th order, in both cases using the last 
50 coefficients. This consistent decrease of the error for increasing m is in 
excellent agreement with our large order predictions. For the case (e) the 
maximum error is 0.15 %. 

The reader may note that A^~\e, 0) is zero for any positive integer e as 
seen in Fig .|To|. Similar zeros in A^~\e,N) are also seen in Fig.[L7|. This is 
due to a new type of supersymmetry for these values of e. We will elaborate 
upon this in the next section. 


6 AT-fold supersymmetry 

In the previous section, we have seen that the Borel singularity of perturba¬ 
tion theory apparently disappears for the first J\f excited states at q = q~ for 
e = A f (A f = 1, 2, 3, • • ■). Hereafter we will call these states isolated states, 
since these states do not have degenerate partners at q = q + (see Fig.[l3|). 
In Table [L. we have listed the perturbative coefficients for several low-lying 
states for J\f = 1. The perturbative correction for the ground state is not only 
Borel summable but in fact identically zero at every order. We also notice 
that the perturbation calculation cannot lift the degeneracy that exists at 


43 





A ] m 



Figure 15: The comparison of the fitting to the perturbative coefficients 
(indicated by dots) and the theoretical prediction (solid line) of 0). 



Figure 16: The comparison of the fitting to the perturbative coefficients 
(indicated by dots) and the theoretical prediction (solid line) of A( + \e, 0). 
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Figure 17: The comparison of the fitting to the perturbative coefficients 
(indicated by dots) and the theoretical prediction (solid line) of (e, N) 


for IV = 1 ~ 3. 
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the zeroth order in g: Ep \l,N,m) and Ep^\l,N — 1 , m) are identical for 
all m. Interestingly similar features prevail in the results of the other values 
of J\f. The results for J\f = 2, 3 show similar characteristics as seen in Table 

|| and |3j. 

The degeneracy reminds us of the non-renormalization theorem of super- 
symmetric theories. This holds that quantum correction at any finite order of 
perturbation theory cannot induce supersymmetry breaking if it is unbroken 
at the tree level, and the present model can actually be cast into a super- 
symmetric quantum mechanics at J\f = 1 [|TJ . The degeneracy between the 
perturbative coefficients, Ep and Ep^l, N — 1, m), is then under¬ 

stood as that between the bosonic and fermionic levels of the perturbatively 
unbroken supersymmetry. In this section we will show that this explanation 
also applies to other A/". We will reveal a new symmetry of the model that 
reduces to the ordinary supersymmetry at Af = 1 and can be understood as 
one of its possible extensions. We dub it A/”-fold supersymmetry. 

Considerations based on this symmetry will further reveal that the per¬ 
turbative series of the isolated levels is convergent with a finite convergence 
radius. The perturbative series of these levels has a much stronger property 
than expected by the lack of the Borel singularity. To begin with, we briefly 
review the supersymmetric quantum mechanics, focusing upon the specific 
realization of the non-renormalization theorem. 


6.1 Supersymmetric quantum mechanics 

We will start by clarifying the supersymmetry that underlies our model at 
Af — 1. To see this, let us define the supercharges |2l|] by 


Q ] = Dif\ Q = 


where 


( 6 . 1 ) 

( 6 . 2 ) 


D—p — iW(q), D Jt =p + iW(q), 
p = —i(d/dq ) is the canonical momentum of q, and W(q) is defined by 

W(q) = q{l- gq). (6.3) 

In the definitions ( |6.1| ), if and if' are the annihilation and creation operators 
for a fermionic degree of freedom. They satisfy = 1 and if^ 2 = if‘ 2 = 

0 . As a result, Q 2 and Q^ 2 vanish: 


Q 2 = Q ]2 = 0. 

The supersymmetric Hamiltonian H is given by 

H = \ {Q*. Q} = \ (p 2 + W 2 (q)) + W'(q) A/>V - \ 


(6.4) 


(6.5) 
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m 

4 4) 

4 } (i),4 +) (o) 

4 } (2),4 +) (l) 

4 } ( 3 ), 4 +) (2) 

4 } (4), 4 +) (3) 

0 

0.5 

1.5 

2.5 

3.5 

4.5 

1 

0 

-3 

-1.2 x 10 1 

-2.7 x 10 1 

-4.8 x 10 1 

2 

0 

-1.95 x 10 1 

-1.41 x 10 2 

-4.665 x 10 2 

-1.098 x 10 3 

3 

0 

-2.70 x 10 2 

-3.330 x 10 3 

-1.5930 x 10 4 

-4.9320 x 10 4 

4 

0 

-5.1791 x 10 3 

-1.0474 x 10 5 

-7.1373 x 10 5 

-2.8885 x 10 6 

5 

0 

-1.2110 x 10 5 

-3.8959 x 10 6 

-3.7318 x 10 7 

-1.9624 x 10 8 

6 

0 

-3.2594 x 10 6 

-1.6258 x 10 8 

-2.1630 x 10 9 

-1.4695 x 10 10 

7 

0 

-9.7888 x 10 7 

-7.3976 x 10 9 

-1.3523 x 10 11 

-1.1804 x 10 12 

8 

0 

-3.2201 x 10 9 

-3.6075 x 10 11 

-8.9708 x 10 12 

-1.0008 x 10 14 

9 

0 

-1.1464 x 10 11 

-1.8643 x 10 13 

-6.2466 x 10 14 

-8.8621 x 10 15 

10 

0 

-4.3816 x 10 12 

-1.0131 x 10 15 

-4.5323 x 10 16 

-8.1376 x 10 17 

11 

0 

-1.7875 x 10 14 

-5.7580 x 10 16 

-3.4085 x 10 18 

-7.7085 x 10 19 

12 

0 

-7.7506 x 10 15 

-3.4088 x 10 18 

-2.6465 x 10 2 ° 

-7.5041 x 10 21 

13 

0 

-3.5604 x 10 17 

-2.0961 x 10 2 ° 

-2.1154 x 10 22 

-7.4855 x 10 23 

14 

0 

-1.7282 x 10 19 

-1.3358 x 10 22 

-1.7367 x 10 24 

-7.6343 x 10 25 

15 

0 

-8.8443 x 10 2 ° 

-8.8090 x 10 23 

-1.4620 x 10 26 

-7.9464 x 10 27 

16 

0 

-4.7629 x 10 22 

-6.0040 x 10 25 

-1.2602 x 10 28 

-8.4300 x 10 29 

17 

0 

-2.6943 x 10 24 

-4.2259 x 10 27 

-1.1113 x 10 3 ° 

-9.1043 x 10 31 

18 

0 

-1.5983 x 10 26 

-3.0695 x 10 29 

-1.0016 x 10 32 

-1.0001 x 10 34 

19 

0 

-9.9264 x 10 27 

-2.2997 x 10 31 

-9.2230 x 10 33 

-1.1167 x 10 36 

20 

0 

-6.4447 x 10 29 

-1.7765 x 10 33 

-8.6724 x 10 35 

-1.2667 x 10 38 


Table 1: The perturbative coefficients of the energy levels E^\l, N,m) de¬ 
fined in Eq. (|5.5|) . In the first row E^\N) is a shorthand denotation for 
E^\l,N,m). 
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m 

4 } (o) 

4 } (!) 

4 } (2),4 +) (o) 

4 } (3),4 +) (l) 

4 } (4),4 +) (2) 

0 

0.5 

1.5 

2.5 

3.5 

4.5 

1 

0 

0 

-6 

-1.8 x 10 1 

-3.6 x 10 1 

2 

0 

0 

-5.1 x 10 1 

-2.55 x 10 2 

-7.14 x 10 2 

3 

0 

0 

-9.09 x 10 2 

-7.227 x 10 3 

-2.7954 x 10 4 

4 

0 

0 

-2.2136 x 10 4 

-2.7102 x 10 5 

-1.4323 x 10 6 

5 

0 

0 

-6.4877 x 10 5 

-1.1942 x 10 7 

-8.5387 x 10 7 

6 

0 

0 

-2.1620 x 10 7 

-5.8676 x 10 8 

-5.6267 x 10 9 

7 

0 

0 

-7.9445 x 10 8 

-3.1262 x 10 10 

-3.9879 x 10 11 

8 

0 

0 

-3.1605 x 10 10 

-1.7757 x 10 12 

-2.9907 x 10 13 

9 

0 

0 

-1.3451 x 10 12 

-1.0635 x 10 14 

-2.3481 x 10 15 

10 

0 

0 

-6.0762 x 10 13 

-6.6658 x 10 15 

-1.9161 x 10 17 

11 

0 

0 

-2.8970 x 10 15 

-4.3483 x 10 17 

-1.6165 x 10 19 

12 

0 

0 

-1.4520 x 10 17 

-2.9405 x 10 19 

-1.4046 x 10 21 

13 

0 

0 

-7.6282 x 10 18 

-2.0554 x 10 21 

-1.2533 x 10 23 

14 

0 

0 

-4.1917 x 10 20 

-1.4817 x 10 23 

-1.1458 x 10 25 

15 

0 

0 

-2.4051 x 10 22 

-1.0998 x 10 25 

-1.0714 x 10 27 

16 

0 

0 

-1.4391 x 10 24 

-8.3939 x 10 26 

-1.0233 x 10 29 

17 

0 

0 

-8.9697 x 10 25 

-6.5820 x 10 28 

-9.9710 x 10 3 ° 

18 

0 

0 

-5.8182 x 10 27 

-5.2988 x 10 3 ° 

-9.9051 x 10 32 

19 

0 

0 

-3.9242 x 10 29 

-4.3774 x 10 32 

-1.0025 x 10 35 

20 

0 

0 

-2.7497 x 10 31 

-3.7096 x 10 34 

-1.0331 x 10 37 


Table 2: The perturbative coefficients of the energy levels (2, iV, rn). 
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4 } (o) 

4 4) 

4 4) 

4 } (3), 4 +) (o) 

4 } (4),4 +) (l) 

0.5 

1.5 

2.5 

3.5 

4.5 

-1 

2 

-1 

-1.0 x 10 1 

-2.5 x 10 1 

1.5 

0 

-1.5 

-1.05 x 10 2 

-4.125 x 10 2 

-4 

8 

-4 

-2.290 x 10 3 

-1.36 x 10 4 

1.3125 x 10 1 

0 

-1.3125 x 10 1 

-6.7676 x 10 4 

-5.9131 x 10 5 

4.8 x 10 1 

9.6 x 10 1 

-4.8 x 10 1 

-2.3888 x 10 6 

-3.009 x 10 7 

1.8769 x 10 2 

0 

-1.8769 x 10 2 

-9.5181 x 10 7 

-1.7008 x 10 9 

-7.68 x 10 2 

1.536 x 10 3 

-7.68 x 10 2 

-4.1525 x 10 9 

-1.0385 x 10 11 

3.2477 x 10 3 

0 

-3.2477 x 10 3 

-1.9476 x 10 11 

-6.7357 x 10 12 

-1.4080 x 10 4 

2.8160 x 10 4 

-1.4080 x 10 4 

-9.7049 x 10 12 

-4.5906 x 10 14 

6.2247 x 10 4 

0 

-6.2247 x 10 4 

-5.0965 x 10 14 

-3.263 x 10 16 

-2.7955 x 10 5 

5.5910 x 10 5 

-2.7955 x 10 5 

-2.8048 x 10 16 

-2.4059 x 10 18 

1.2718 x 10 6 

0 

-1.2718 x 10 6 

-1.6112 x 10 18 

-1.8330 x 10 2 ° 

-5.8491 x 10 6 

1.1698 x 10 7 

-5.8491 x 10 6 

-9.6315 x 10 19 

-1.4386 x 10 22 

2.7148 x 10 7 

0 

-2.7148 x 10 7 

-5.9791 x 10 21 

-1.1606 x 10 24 

-1.2701 x 10 8 

2.5402 x 10 8 

-1.2701 x 10 8 

-3.8483 x 10 23 

-9.6071 x 10 25 

5.9828 x 10 8 

0 

-5.9828 x 10 8 

-2.5649 x 10 25 

-8.1491 x 10 27 

-2.8353 x 10 9 

5.6706 x 10 9 

-2.8353 x 10 9 

-1.7687 x 10 27 

-7.0762 x 10 29 

1.3508 x 10 10 

0 

-1.3508 x 10 10 

-1.2610 x 10 29 

-6.2853 x 10 31 

-6.4664 x 10 10 

1.2933 x 10 11 

-6.4664 x 10 10 

-9.2892 x 10 3 ° 

-5.7075 x 10 33 

3.1087 x 10 11 

0 

-3.1087 x 10 11 

-7.0678 x 10 32 

-5.2964 x 10 35 


Tabic 3: The perturbative coefficients of the energy levels E^\s, N,m) 
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All the states are labeled by the quantum number of the fermion number 
operator Nf — ip' ^The operator Qi transforms a bosonic state |b) (with 
Nf — 0) to a fermionic one |f) (with Nf = 1), and Q does vice versa. By 
construction, Q and Q t commute with H 


[H,Q] 


H.gt 


0 . 


( 6 . 6 ) 


Therefore, the Hamiltonian is invariant under the transformations generated 
by the supercharges. 

The relation between this supersymmetric Hamiltonian and our Hamil¬ 
tonian becomes obvious if we introduce a matrix representation of 



(6.7) 


In this representation, the bosonic state |b) and the fermionic one |f) are 
written as 


|b> = 


0 

vp(-) 



( 6 . 8 ) 


and the supersymmetric Hamiltonian is decomposed into two components: 


H = 




(6.9) 


Here we denote the fermionic component by H + and the bosonic by H_ which 
are given by 

H- = H + = l -DD\ (6.10) 

or more explicitly 

(6 - n) 

One can readily check H + — 1/2 is indeed the Hamiltonian of our model 
with J\f = 1, i.e., the one derived from the action 0-@) with e — 1. 
The Hamiltonian H_ — 1/2 becomes H + — 1/2 as a result of the mirror 
transformation with respect to the q = 1/2 g axis. Therefore, the dynamics 
of our model with J\f = 1 are same as that of the supersymmetric Hamiltonian 
H. 

In considering the eigenvalues of H, we will work on the decomposed form 
of the equation, i.e., 

# ± T (±) = £WT (±) . (6.12) 

The perturbative energy levels E calculated in the previous section were 
those of the well-local states in the left well of H + — 1/2 and of the states 


50 




A 




H + 



Figure 18: The correspondence between the states in the left and right wells 
of the potential (|2.2|) with e = 1 and the states in the left wells of H± — 1/2. 


in the right well. It is convenient to identify the latter states with the well- 
local ones in the left well at q ~ 0 of iF_ — 1/2 (see Fig.[L8]). These states 
correspond to “bosonic” states, while the well-local states in the left well of 
H + — 1/2 correspond to “fermionic” states. 

At g — 0, H reduces to that of the supersymmetric harmonic oscillator. 
The eigenvalues E^> and eigenfunctions can p e solved immediately. 
They are labeled by an integer N(= 0,1, 2,...) representing the level of exci¬ 
tation. The solution for is 

^\q) = H N {q)e-« 2 / 2 (6.13) 

by the use of the N-th Hermite polynomial H N (q ) = (— l) N e q ~ -j^e~ q2 ; the 
corresponding eigenvalues are * = N and — N + 1. The ground 
state Tq ^ is of zero-point energy and the supersymmetry is not broken. The 
degeneracy between the fermionic and bosonic spectra due to the unbroken 
supersymmetry is manifest in the pairs of E and E^. 

Now let us examine a case where g is a nonzero quantity. When one 
applies perturbation theory to Eq. (|6.12|) by expanding the wave functions 
and the eigenvalues in g, one gets 

1 OO 

4? = -F + E E^\l,N,m)g 2m , (6.14) 

Z m =0 
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where E\^ (1, N,m) is exactly what we have calculated in the previous sec¬ 
tion, Eq.( |5.5|) . If we recall the identification of the state in the right well of 
H A) and the one in the left of H^~\ it will be seen that the perturbative 
coefficients of E jy are given by those of Ep ■’(l, N) calculated in the previous 
section. The disappearance of the correction to Eq ^ and the persistence of 
the degeneracy between and E^ are now seen in the context of the 

unbroken supersymmetry of H in perturbation. 

A direct demonstration of the non-renormalization theorem can be ob¬ 
tained for the present case, by considering the function obtained from -DTg = 
°, 

'M?) = exp £ dq'W(q')^ = exp j • (6-15) 


Since this function starts with v h[ ) ~' ) (g) in the expansion in g and obviously 
satisfies H_^/q = 0, it mimics the zero-energy state of H_, whose existence 
is the indicator of unbroken supersymmetry. 

The function Tq is not normalizable, and cannot, therefore, be accepted 
as the physical zero-energy state. Thus, the supersymmetry is dynamically 
broken at any non zero values of g [|2T|. The perturbation theory, however, 
cannot detect it, for the unnormalizability does not appear at any finite order 
in g in the expansion of To• (There can be no zero-energy state in H + either, 
for the spectrum of H + is necessarily the same as that of i/_.) 

Since the perturbation cannot break the supersymmetry, the degeneracy 
of the excitations must persist in the perturbation. Thanks to the simplicity 
of the quantum mechanics, we can check this directly. The explicit contents 
of the supersymmetry relation (3.6) are the two equations 


DH- = H + D , (6.16) 

D ] H + = H_D ] . (6.17) 


We will focus on Eq. (|6.16| ) in the following discussion. This holds that if 
tf(-> is an eigenfunction of H_ with the eigenvalue E 1 '- 1 then ThU ■ 1 is that 
of H + with the same eigenvalue. The operators D and D' are written as 

D = (-i)e (OT3/3) ae“ (OT3/3) , D ] = ie^" 3 ^ a) e^ q3/3 \ (6.18) 


where 


a = 


d_ 

dq 


(6.19) 


is the annihilation operator in the harmonic oscillator up to an unimportant 
factor \/2, [a, ad] = 2. The operator D has a as the leading term in its 
expansion in g\ D — a + 0(g). Thus if the state is for the (N + 1)- 
th excitation and starts with H N+1 (q)e~ q / 2 in the perturbative expansion, 
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D'i’E) starts with H]y(q)e ~ q2 / 2 and represents the iV-th excitation of H + 
This observation suggests the degeneracy in the perturbative energy levels. 
To be more specific we can rewrite the two Hamiltonians as 


H + = (l/2)e (9 ' ?3/3) a(a t -2^ 2 )e- (9 ' ?3/3) , 
H_ = (l/2)e i9q3/3) (a^ -2gq 2 )ae 


—( 99 3 / 3 ) 


( 6 . 20 ) 

( 6 . 21 ) 


In conducting the perturbative analysis of we need not necessarily 

expand in g directly. Since the asymptotic expansion is unique, we 

may adopt another convenient expansion scheme for the eigenfunction as 
long as it is well-defined. Looking at Eqs. Q6.18| ), (|6.20|) , and ( |6.21|) we can 
expand = e _ ( £ rt 3 / 3 )\l/ ( ~) j n g ; 


LXJ 

= H N+l (q)e-'! 2 + V g">«> lm) 


N+l 


(<?)• 


( 6 . 22 ) 


m= 1 


The eigenvalue equation for now reads 


2 ( a ^ 2gq 2 ^) a'I'jvq-i — -^ir+i^jv+n 

By multiplying a on both sides of Eq. (|6.23 ), we obtain the following: 


(6.23) 


-a (a f - 2 gq 2 ^j 


a<h 


(-) ' 
N+l 


— rj N +1 


a<f> 


(-) 
N +1 


(6.24) 


The function [a^^+i] starts with H N (q)e ~^ q2 / 2 - ) when expanded in g. 
Eq. (|6.24j) shows that = e (99 3 / 3 )[ a $w)p as ^j ie same perturbative series 

as Ej^li to the eigenvalue in H + . This proves that the (N + l)-th excitation 
of H _ and the iV-th of H + have the same perturbative series. 

We can now turn to the other values of A/". The first 20 coefficients of 
E^' 1 for J\f = 2 are listed in Table and those for J\f = 3 are listed in Table 
|3|. The vanishing coefficients at J\f = 2 for the two isolated states draw our 
attention first |] The perturbative degeneracy remains, but this time between 
Ep ^ (A/”, N + A/", m) and Ep + ^ (A f, N, m ). In the next section, we will present 
the general proofs of these properties of the perturbative coefficients. 


6.2 Non-renormalization theorem 

We will start with the degeneracy, and define a pair of two Hamiltonians 
H±(J\f) with Af as a parameter, 

H ±W = + i 9 2 (l - gq? T V (gq - 5) . (6.25) 

6 Supersymmetry provides an explanation for the vanishing coefficients for only the first 
excited state [p~8|. 
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The Hamiltonians H±(Af) are related to the action ( | 2 . 1 | )-( | 2 . 2 | ) again with 
e = Af. We denote the wave function and the eigenvalue of the iV-th excited 
state in the left well of H±(Af) by and E^\ as in the case of the 
previous subsection. The degeneracy between the N + Af-th excited level at 
the right well and the iV-th at the left in Table |2] and |3] is realized between 
the eigenvalues of the left wells of H±(J\f). The Hamiltonians are also 
rewritten as 

H+{U) = \DD' (6.26) 

H-(K) = \d'D + (JV - 1 ) ( 9 ? - ±) . (6.27) 


As in the supersymmetric case, these two Hamiltonians transform to each 
other by mirror transformation and thus describe the same quantum me¬ 
chanics. 

The point of the explanation concerning the degeneracy in Af = 1 was 
Eq. (|6.16|) , and the fact that the operator D reduces the level of excitation by 
1. What we have now is degeneracy between E jv+jv' an d M v"^ Correspond¬ 
ingly, we find H±(Af) satisfies 


D*H-{. AT) = H + (Af)D Ar , (6.28) 

D^H + (Af) = H_(AT)D ]M . (6.29) 


The proof is as follows. We first note that the equation 

H + (Af)D = DH + (Af-2) - i(Af - l)g, (6.30) 


is proved from Eqs. (|6.2|) and ( |6.26|) . We examine Eq. (|6.28|) . The other will 
be proved in a similar way. We start with H + (Af)D^ and move each D to 
the left using Eq. (|6.30| ). After moving all D 7 s we obtain 


M-l 

H + (Af)D M = D m H + (AT - 2Af) - igD^- 1 (Af - 1 - 2k). (6.31) 

k =0 


The terms in the summation over k on the right-hand side cancel each other 


and H + {—Af) = H_(J\f). This proves Eq. Q6.28Q . 

Once we have obtained Eq. ( |6.28|) , the perturbative degeneracy is proved 
in a similar manner to the ordinary supersymmetric case. In terms of a and 


a f , 

H ± {Af) = e (9q3/3) h ± e~^ q3/3) , 


(6.32) 

where 

h+ = (l/ 2 )a(a f - 2 gq 2 ) - (Af - 1) ( 

[" ~ s) ’ 

(6.33) 


h _ = (l/ 2 )(a t — 2gq 2 )a + (Af — 1) ( 


(6.34) 


54 
















Substituting Eqs. (|6.18| ) and ( |6.32| ) into Eq. 

a^h- = h + a^ 


leads 


(6.35) 


This establishes that the perturbative series of the (N + A/")-th level in Il¬ 
ls the same as that of the iV-tli level in h + . In turn it establishes the per¬ 
turbative degeneracy between E and The wave function of 

the (N + A/")-th state of H^~\,AT) is related to the wave function of the 
iV-th excited state of (JV): 


D n Tl7 ) 

^ ^N+M *N ■ 


On the contrary, the isolated states should satisfy 

= 0 (./V = 0,1,2, • • •, A/” — 1), 


(6.36) 


(6.37) 


which is the analogue of the ground state condition of the supersymmetric 
case. Since is rewritten as 


l-)M — g-(9 2 / 2 )+(9<? 3 /3) 


. d ' 

_ 

dq / 


M 


3 (q 2 /2)-(sQ 3 /3) 


the solutions for Eq. (|6.37|) have the form 




P = /(«)« p 


(6.38) 


(6.39) 


where f(q ) is a polynomial of q of order A/" — 1 at most. Again the unnormal 
izability of ; disappears once it is expanded in g, as in the case of Af = 1 
Plugging Eq. ( |6.39| ) into 


N i 


(6.40) 

(6.41) 

(6.42) 


we obtain 

f" + 2 (gq 2 - q)f + [—2< 7 (A'” - l)q + 2E + N - 1 ] / = 0. 

Let us write the polynomial / explicitly as 

M-l 

Hq) = a ^ n - 

n =0 

Only polynomials of the order of (A/"— 1) at most can be a consistent solution 
for ( |6.41| ): Otherwise the terms 2 g[q 2 f — (A f — 1 )qf] in ( |6.41| ) forces f(q) to 
have much higher powers.^] 

7 We note that the equation (6.41) is rewritten in the form [— (1/2) (J - ) 2 — gJ + + J°]f = 
Ef where J~ = d/dq , J° = q{d/dq) — (1/2)(A/"— 1), and J + = q 2 (d/dq) — (Af — l)q are the 
SL(2,Q) generators for the Af-dimensional representation. This is the eigenvalue problem 
equation of the quantum top, known in connection with the quasi-exactly-solvable problem 
j22j. This also accounts for the fact that the solution / is obtained in polynomials of the 
order Af — 1. 
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It is convenient to write the resulting equation for a n in a matrix form 


M„{E) • a = 0, (6.43) 

where a T = (ao,ai, The (i,j )-element of is given by 

= 2g(i — Af)§ij + i + (2E-\-J\f — 1 — 2i)hy + (i + 2)(i + l)<5jj_ 2 , (6.44) 
where i and j run from 0 to A/" — 1. The matrix M^/'(£') has the explicit form 
M W{E) = 

( 2E+N-l 0 2 0 ^ 

2g(l - J\f) 2E + J\f — 3 0 6 

0 2g(2 — M) 2E+M-3 0 

0 0 2g(3- Af) 2E + M-7 

V ; / 

(6.45) 

at the upper left corner and 

M n {E) = 


V 


2E — AT + 7 0 (A/” — 2) (A/” — 3) 0 

— 6 g 2E-A7 + 5 0 (A/” - 1)(A/'- 2 ) 

0 -4 g 2E-U + 3 0 

0 0 -2 g 2E-N+1 


at the lower right. 

The condition that there exists a non-trivial solution for a is 


(6.46) 


det M n {E) = 0. (6.47) 

Writing this condition explicitly for several lowest values of A/", we obtain 


A/" = 1 : 
A7 = 2 : 
M = 3 : 
A/" = 4 : 
A/” = 5 : 


E = 0, 

E- - 
2 


e + - 2 )=o, 


E(E-l) (E+l) + 2g 2 = 0, 


E 


E — 


E + 


E + -)+ 12 Eg 2 = 0 , 


E[E - 2)(E - 1 )(E + 1 )(E + 2) + 6(7 E 2 - 52 )g 2 = 0. 


(6.48) 
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These algebraic equations for E of the order J\f determine the perturbative 
eigenvalues of the J\f isolated states. The solutions are consistent with those 
we expect as the energy levels in perturbation theory. They reproduce the 
energy levels of the first A f lowest states of the harmonic oscillator at g — 0, 
which are E = -(JV/2) + (1/2), -(A/"/2) + (3/2),..., (A/"/2) - (1/2). The 
perturbative corrections are identically zero at any g 2 for J\f = 2 as well as 
J\f — 1. We have also checked that the equation at other values of J\f (A f > 3) 
correctly regenerates the perturbative coefficients obtained in the previous 
section once it is solved in powers of g 2 . 

The very existence of the algebraic equation (|6.47|) for E tells us that the 
perturbative series for the isolated energy levels are convergent, not simply 
Bor el summable. 

Finally, we will combine the relations ( |6.28| ) and (|6.29 ) into a compact 
form, which suggests the name “A/"-fold supersymmetry”, and will present an 
interesting relation between the Hamiltonians H±(AT) and a new Hamiltonian 
that will be introduced naturally as a direct extension of the supersymmetric 
Hamiltonian. 

We will define the A/"-fold supercharge in 2 x 2 matrix notation by 



Qat ~ 


Qat 


0 0 

{D^) N 0 


(6.49) 


(6.50) 


These obviously satisfy 

Qm = QV = 

and induce transformations between the “bosonic” state | B) and the 
“fermionic” state \E)\ 


\B) = 


0 


l^) = 


g>(+) 


x _ ^ 0 
The relations ( |6.28|) and (|6.29|) now become compact as 


[Hat, Qat ] = -Hjy, Q\r 


= 0 


where H^r is defined by 


H M = 


H + (JV) 0 
0 H_( AT) 


in terms of the Hamiltonians H±(J\f) in Eq. 

Hamiltonian is constructed from the supercharges Q and Q': 

D n D^ m 0 
0 D^D U 


Km=\ {Qn, Q ] Af] = ^ 


(6.51) 


(6.52) 


(6.53) 


The supersymmetric 


(6.54) 
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which we call the A/"-fold supersymmetric Hamiltonian. The Hamiltonian 
Hy commutes with Qy and Qy. We can define the states 

|B 0 > = ^h), n = (6.55) 

using the perturbative isolated states of H_ (J\f). These states satisfy 

Hy|£ 0 > = 0, (6.56) 

and play the role of the A F vacua of Hy. The symmetry generated by Qy 
and Q^- is A/"-fold of the ordinary supersymmetry in the sense of the number 
of perturbative vacua and the change in the level of excitation. 

By the definition of Hy and Eq. (|6.52| ), it is obvious that 

[Hy, Hy] = 0. (6.57) 

Therefore Hy is a function of Hy. By counting the highest orders of d/dq 
in both operators, we find that Hy must be a polynomial of Hy of order 
A f. Since the perturbative isolated states satisfy Eqs. (|6.47| ) and (|6.56|) , we 
conjecture that 

Hy = ^ det My(Hy). (6.58) 

In the above, we have fixed the overall constant by comparing the highest 
powers of d/dq on both sides. We have checked this conjecture explicitly by 
hand to the value Af = 3, and algebraically on a computer to A f — 10. 

The non-renormalization theorem in the A/"-fold supersymmetric quan¬ 
tum mechanics guarantees the absence of any perturbative correction to the 
energy of the ground states in Hy and the persistence of the degeneracy in 
perturbation theory. 


7 Numerical verification of the prediction of 
the energy spectrum 

In this section we will numerically verify the prediction with regard to the 
energy spectrum made from the valley method. The result (|4.16|) obtained 
in Section 4 is for the case of no perturbative degeneracy, or even for the 
isolated states in the presence of degeneracy. The result ( |4.18| ) is for the 
perturbatively degenerate states. In principle we can use these to check the 
validity of the valley method by a comparison with exact numerical energy 
levels. There is, however, a great difficulty for general values of e. The valley 
result represents only the non-perturbative part and needs to be combined 
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with the perturbative predictions. This is important for numerical studies, 
since the perturbative part has in general a much larger numerical value 
than does the non-perturbative. The problem is that we know only the 
leading Borel singularity of the perturbation series. After canceling it with 
the imaginary part of the non-perturbative contribution we obtained, we are 
still left with non-leading Borel singularities in the perturbation theory. This 
prohibits us from making numerical predictions for the energy levels. 

The model with e = A f has a great advantage in this respect. For the 
isolated states we can obtain the perturbative contribution exactly by solving 
Eq. (|6.47|) . Remember that the equation generates all the perturbative coeffi¬ 
cients and, thus, it contains complete information on the perturbative contri¬ 
bution. Correspondingly, the valley prediction for these states derived from 


Eq.( 4.16 ) is free from any imaginary element. Regarding the perturbatively 
degenerate states, their energy splittings have no perturbative corrections 
because of the A/”-fold supersymmetry. 

In this section we will carry out the numerical calculation for the energy 
spectra in the case of Af = 0 , 1 , 2 , 3, and compare them with the valley pre¬ 
dictions. We label its energy levels by a single integer k , E k (k = 0,1, 2,...), 
counting from bottom. That is, the energy of the ground state, the first 
excited state, and so on, are denoted by E 0 , Ei, ..., respectively. They thus 
include the isolated states localized at the lower right well. We will first make 
the correspondence with the result in the previous sections and identify the 
quantities on which we will perform the comparison. 

For A f isolated states, k — 0,1 ,...,Af — 1, we use E^' 1 for denoting the 
perturbative contribution to the energy, which is the solution of Eq. ( |b.47| ). 
In fact they are zero at J\f = 1 and just —1/2 or 1/2 at J\f = 2. In the case 
of Af = 3, the equation reads 

E(E — 1)(E + 1) + 2g 2 = 0. (7.1) 

This equation has three real roots up to g 2 = \/3/9 and each of them repre- 
sents E k . We will first solve the equation numerically at various values of 
g, and then we will define the non-perturbative correction by 


(P) 


AE k = E k — E k 
We can compare this with the valley prediction 


A 4 V) = a 2 - 


- 1 ) 


fc+i 


—Af+2k 


k\ 


-T(Af - k ) 


(7.2) 


(7.3) 


which is derived from Eq. (|4.16|) by the insertions e = Af and IV_ = k. 

For the perturbatively degenerate energy levels, k = A/", Af + 1,..., we will 
consider the energy splitting defined by the following equation: 


A E 


k fc+1 


= E 


fc+i 


E k . 


(7.4) 
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8 

Figure 19: The comparison in absolute values on the logarithmic scale be¬ 
tween the instanton prediction and the exact numerical result for the energy 
splitting at J\f = 0. The solid lines stand for the instanton predictions, while 
the dots show the corresponding exact numerical results. 


The valley prediction in relation to this splitting is 


A-Efhi = 2a. 


\ p-A0/2]!p+A0/2]! W 


(7.5) 


as derived from Eq. (|4.18|) by the insertion of N± = (k =|= A/")/2. 

The result of the comparison is plotted in rigs. [TH| - |T2| . The valley predic¬ 
tion monotonically approaches to the exact values as g decreases for all levels. 
This fact clearly shows that the results (|4.16|) and ( |4.18|) obtained on the ba¬ 
sis of the valley method are indeed the leading terms of the non-perturbative 
correction to the energy spectrum. 

There are various sources of corrections to the leading term that could 
account for the discrepancy in the figures. First is the fact that the potential 
(|2.2|) ceases to have two minima at g larger than \J v^3/(18A/"), that are g > 
0.31 at Af = 1, 0.22 at 2, and 0.18 at 3, respectively. For those values of g, the 
valley-instanton is no longer a solution of the valley equation and the valley 
prediction becomes unreliable. Thus it is natural for the valley prediction 
to have a relatively large discrepancy at these values as above. Even for the 
smaller g , there must be higher-order corrections from fluctuations in the 
valley background. The inclusion of those corrections will generate an 0(g 2 ) 
correction relative to the leading term. 

In Tables £|-|(|, we list the specific values of A Ej., AEkk+i, A E\ 
and the relative discrepancy 5 between them, 


(V) A n(V) 

A^fefe + 1 , 


5 = 


A£( y ) - A E 


A E 


(7.6) 
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Figure 20: The comparison between the valley prediction and the exact nu¬ 
merical result for the non-perturbative corrections to the isolated energy level 
and splitting at A f = 1. The solid lines stand for the valley prediction, while 
the dots show the corresponding exact numerical result. 



g 

Figure 21: The comparison between the valley prediction and the exact nu¬ 
merical result at A f — 2. We use dashed lines to represent the valley pre¬ 
diction at g larger than 0.22, as the valley prediction has no justification for 
those values, as explained in the text. 
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Figure 22: The comparison between the valley prediction and the exact nu¬ 
merical result at AT = 3. We use dashed lines to represent the valley pre¬ 
diction at g larger than 0.18, as the valley prediction has no justification for 
those values, as explained in the text. The critical value of g for the equation 
(17.1|) to cease to have real solutions is well beyond the scope of this figure. 


Because of the accumulation of round-off error in the calculation of the en¬ 
ergy, it is difficult to carry out the numerical calculation for smaller g than 
shown in the figures and tables. The discrepancy is not necessarily small, 
even at the lowest g. Nevertheless the discrepancy decreases as g decreases 
and the behavior of the decrease seems consistent with the expected depen¬ 
dence on g 2 . The legitimacy of our valley prediction as the leading term of 
the non-perturbative contribution is solid. 

The excellence of the agreement of the valley prediction and the exact 
result for the ground state energy in both AT = 1 and 2 is very striking. 
This may indicate a new sort of non-renormalization theorem, even in the 
non-perturbative sector. 

Finally, we note that these non-zero values of A E are the quantitative 
measure of the dynamical breaking of the AAfold supersymmetry. 


8 Discussions 

1. First we would like to outline the results of this paper. We have exam¬ 
ined the geometrical structure of the configuration space of the path 
integral, and shown that in one-dimensional quantum mechanics: 

(a) There are valleys in the configuration space which connect the 
vacuum and the instanton-like configuration. (Section ^J) 
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Table 4: The values of non-perturbative part of the energy spectrum at A f = 
1 ; A E is from the exact numerical calculation, Athe valley prediction, 
and <5 is their relative discrepancy S = |(A— AE)/AE\. 
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g 

AEq 

AEq V) 

h 

0.14 

-6.438xl0“ n 

-6.413x10“®® 

0.00387 

0.16 

-4.530x 10“ 9 

-4.510xl0“ 9 

0.00445 

0.18 

-8.825x10-® 

-8.776x10“® 

0.00569 

0.20 

-7.706xl0- 7 

-7.651xl0“ 7 

0.00713 

0.22 

-3.967xl0 -6 

-3.932xl0“ 6 

0.00877 

0.24 

-1.421xl0 -5 

-1.406xl0“ 5 

0.0106 

0.26 

-3.933x 10 -5 

-3.884xl0“ 5 

0.0126 

0.28 

-9.018x 10 -5 

-8.884xl0“ 5 

0.0148 

0.30 

-1.795xl0 -4 

-1.764xl0“ 4 

0.0170 

g 

AEi 

A E[ V) 

h 

0.12 

1.788xl0- y 

1.956xl0“ y 

0.0942 

0.14 

5.883xl0 -7 

6.678xl0“ 7 

0.135 

0.16 

2.317xl0- 5 

2.753xl0“ 5 

0.188 

0.18 

2.659X10" 4 

3.344xl0“ 4 

0.258 

0.20 

1.423x10-® 

1.913x10“® 

0.344 

0.22 

4.658x10-® 

6.714x10“® 

0.441 

g 

AE 23 

A W 

h 

0.08 

1.260x10“® 

1.527x10“® 

0.212 

0.10 

6.689xl0“ 5 

9.220xl0“ 5 

0.378 

0.12 

5.211x10“® 

8.687x10“® 

0.667 

0.14 

5.139xl0“ 2 

1.179x10“® 

1.29 

0.16 

1.576X10- 1 

5.797x10“® 

2.68 

0.18 

2.765X10- 1 

1.596 

4.77 

0.20 

3.826x10“® 

3.093 

7.08 


Table 5: The values of non-perturbative part of the energy spectrum at 
A/" = 2. Since the numerical agreement of Awith the exact value is 
excellent, we think it permissible to write them up to g — 0.3, although the 
valley-instanton ceases to be a solution of the valley equation at g > 0.2. 
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g 

AEq 

A E^ 

5 

0.15 

1.305X10" 11 

1.483x 10 -11 

0.137 

0.16 

1.002xl0” 10 

1.155xlO” 10 

0.152 

0.17 

5.553xlO" 10 

6.508xl0~ 10 

0.172 

0.18 

2.383X10' 9 

2.843xl0~ 9 

0.193 

g 

AEi 

A E[ V) 

5 

0.13 

-4.738xlO~ iU 

-4.324xlO” iU 

0.0874 

0.14 

-7.289xl0 -9 

-6.544x 10 -9 

0.102 

0.15 

-6.639xl0 -8 

-5.859x 10 -8 

0.117 

0.16 

-4.068xl0 -7 

-3.523x 10 -7 

0.134 

0.17 

-1.836xl0 -6 

-1.558x 10 -6 

0.151 

0.18 

-6.527xl0 -6 

-5.417x 10 -6 

0.170 

g 

AE 2 

A W> 

5 

0.11 

1.868X10' 9 

2.362xl0~ 9 

0.265 

0.12 

1.019xl0“ 7 

1.358xl0 -7 

0.333 

0.13 

2.144xl0 -6 

3.028xl0 -6 

0.412 

0.14 

2.257xl0~ 5 

3.407xl0' 5 

0.510 

0.15 

1.420xl0“ 4 

2.315xl0 -4 

0.630 

0.16 

6.040 xl0~ 4 

1.075xl0~ 3 

0.780 

0.17 

1.900xl0~ 3 

3.731xl0~ 3 

0.964 

0.18 

4.732xl0~ 3 

1.032xl0~ 2 

1.18 

g 

AE 3 4 

A W 

5 

0.08 

1.140 x10 _y 

1.559xlO~ Y 

0.367 

0.09 

1.524xl0~ 5 

2.301xl0 -5 

0.509 

0.10 

4.415xl0~ 4 

7.528xl0~ 4 

0.705 

0.11 

4.657xl0“ 3 

9.275xl0~ 3 

0.991 

0.12 

2.402X10' 2 

5.910xl0~ 2 

1.46 

0.13 

7.170xl0 -2 

2.378X10" 1 

2.32 

0.14 

0.143X10” 1 

6.878X10" 1 

3.81 


Table 6: The values of non-perturbative part of the energy spectrum at 
N = 3. 
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(b) The perturbative and non-perturbative contributions in the path 
integral can be separated in terms of the valley. (Section |3|) 

(c) In the separation, the imaginary part of the non-perturbative con¬ 
tribution and the Borel singularity of the perturbative contribu¬ 
tion cancels trivially. (Section |0) 

(d) Bogomolny’s trick is proved to be sound, thus enabling us to go 
beyond the dilute-gas approximation. (Section and Section 

H) 

(e) Based on the above separation, the large order behavior of the 
perturbation series is predicted. The prediction reproduces the 
perturbation theory excellently. (Section |T3| and Section ||) 

(f) A new type of supersymmetry in quantum mechanics is found. 
Non-renormalization theorems of the A/"-fold supersymmetry hold. 
(Section |(|) 

(g) The energy spectrum calculated by our valley calculus reproduces 
those derived directly from the Schrodinger equation. (Section [?]) 


2. The asymmetric double-well potential was also treated by Balitsky and 
Their calculation, however, utilized a complex solution, 


Yung [23 


which requires the deformation of the path-integral in the infinite di¬ 
mensional configuration space. This procedure is difficult to justify, 
and furthermore it obscures the analytic continuation in the complex 
g 2 -plane. On the other hand, our analysis involves only the real con¬ 
figurations, and the analytic continuation is well-defined in the finite 
dimensional collective coordinate space. 


3. To parameterize the valley, we have used R, which corresponds to the 
relative distance between / and /. But as was denoted in the footnote 
to Section [|, the definition of R is not simple. The simplest variable to 
parameterize the valley is the eigenvalue A in the valley equation (|2.5|) . 
In further research on valleys, one might use A rather than R. 


4. Our analyses have been done in the topologically trivial sector, but the 
extension to the topologically non-trivial sector is straightforward. For 
example, we can construct a valley in the one instanton sector which 
consists of the n + 1 instantons and the n anti-instantons. The valley 
equation is the same as in the topological trivial sector, and the only 
change is the boundary condition of q at r = ±oo. 

5. In Section |3]3], we have predicted the large order behavior of the per¬ 
turbative series for the square of the absolute value of the normalized 
wave function at q = q± (see Eq. (|3.38|) ). We have not tested this pre¬ 
diction, but we note that the prediction is consistent with the result in 
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Section |j. In Section ||. we have shown that if e = AT the perturbative 
ground state (or more generally first J\f perturbative excited states) at 
q — q_ can be constructed exactly. The state is normalizable at any 
finite order of the perturbation of g, but if we sum up the contributions 
at all orders, the state becomes minormalizable. This unnormalizabil- 
ity indicates the non-Borel-summability of the the perturbative series 
for the square of the absolute value of the normalized wave function at 
q = q_. 


6 . In Section |7j. we have found in Figs.|(J-|2| that for the first (A f — l)-th 
energy levels our valley calculus agrees excellently with the numeri¬ 
cal calculation. In spite of the fact that we have only taken into ac¬ 
count the lowest order corrections in the background of the valley, and 
neglected the higher order corrections, the agreement is better than 
expected. This may suggest another non-renormalization theorem for 
the higher order corrections on the valley calculus. Indeed, for the 
instanton calculus in the supersymmetric QCD, similar types of the 
non-renormalization theorem are known to hold |24 . 


7. The techniques of this paper may be applied to the model in quantum 
held theory. In quantum held theory, there are three types of singu¬ 
larities in the Borel plane, which are known as the IR renormalon, UV 
renormalon, and the instanton singularities [ffl, [15|. The former two 


singularities are outside the scope of this paper, but the instanton sin¬ 
gularity is manageable by a straightforward extension of our method. 
A promising application of our method is to processes with huge num¬ 
bers of particles. The number of Feynman diagrams of these processes 
is large and the instanton type singularity is closely related with this 
fact [^5, 


15 


8 . One of the interesting processes with high multiplicities is the elec- 
troweak baryon and lepton number violating process at high energy 
scattering |2(| |27|. The process is the tunneling effect, and the height 
of the tunneling barrier is 0(10) TeV. The amplitude of the process 
at that energy is naively expected to be enhanced, and the process in¬ 
evitably contains a huge number of particles (massive gauge and Higgs 
bosons) [29|, [k| . A plausible calculation does not yet exist. 

Based on the optical theorem, Porrati and later Khoze and Ringwald 
proposed a method to treat a huge number of particles |3T], [32|. The 
problem with their method is that the optical theorem itself does not 
distinguish the baryon number violating process (tunneling effect) and 
the baryon number preserving one (the perturbative effect). The exten¬ 
sion of our method to the electroweak theory might make it possible to 
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separate them, although the detailed analysis of the geometrical struc¬ 
ture of the theory is needed. The valley instanton in the electroweak 
theory has already been constructed in Ref. |33 . 


More detailed arguments about this subject are given in Ref. [34 


9. The extension of our method to string theory is an interesting possi¬ 
bility. Unfortunately, it is difficult at present because of the lack of a 
satisfactory non-perturbative formulation of string theory (e.g. string 
held theory). But the consideration of the large order behavior of string 
theory [|35| and the existence of the solitonic object suggests that a 
geometrical structure similar to that of our model exists in the config¬ 
uration space of the non-perturbative formulation of the string theory. 
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A Some notes on the valley method 


In this appendix we will note some of the important properties of the valley 
method, which have not been explicitly described in previous references. 

Firstly, we note its relation to the ordinary collective coordinate method. 
In our context, this corresponds to the cases when there are non-trivial so¬ 
lutions of equations of motion. That is, cases where F(r) = 0. These cases 
can be treated by a limiting procedure and the valley method reduces to the 
usual collective coordinate method. This is done in the following manner: 
When F(t) = 0, both the numerator and the denominator of the Jacobian 
A[ip a \ of Eg. ( |2.12|) are zero. This quantity can be obtained by introducing 
a small perturbation to the action that renders F(t) A 0 and making the 
perturbation vanish at the end. Let us take S[q] —> S[q] + pW[q\. Since 
the eigenvalue A is zero p 7 ^ 0, A is of order p and the valley solutions are 
expanded as q = qo + pq\ + • • •, F — pFi -|— ■. The p-expansion of Eq. (|2.17p 
and Eq. (|2.18|) yields the following: 


/ 


dr' 


t 2 S[q 0 } 

5q 0 (T)6q 0 (T') 


■5%o] _ fl 

Sqo(r) 

(A.l) 

-Fi{t') = 0. 

(A.2) 


The first equation is the original equation of motion. The second gives 
F\ (X qo, which proportional constant is to be determined at the higher-order 
equations. The gradient vector G(t) is independent from this proportional 
constant: 


G(t) = 


FAt) 


? o ( t ) 


(A.3) 




By choosing the valley parameter a as the translation coordinate of q , o( r ) ! 
the Jacobian is found to be the following: 


A[q a ] = J qo(T)G(r)dr = J qfidr'. 


(A.4) 


This is the Jacobian of the ordinary collective coordinate method & 

The second is the extension to the multi-dimensional valley, which will 
be needed when there are multiple dangerous eigenvalues. The definition of 
the .D-dimensional valley is as follows: 


<*%] 

Sq(r) 


D 


dr 


5q(r)5q(T' 


1=1 

(A.5) 

A iFi(r) (i = 1, •••,£>), 

(A.6) 
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The Gaussian integration is performed in the subspace defined by 
/ (f) (r)Gj(r )dr = 0 for i = 1 |7|], where G 3 is the normalized func¬ 

tion of Fi. The Jacobian of the collective coordinate volume element Ylidcii 
is given by the determinant of 

Cij = I dr'^G-G j (r% (A.7) 

where the valley is parametrized by the D collective coordinates cq. 


B The determinant of the valley-instanton 


In this appendix we will calculate the determinant for a valley-instanton, 
utilizing the fact that it has a zero eigenvalue, as is proven in Section 2.2. 
We will consider the valley-instanton located at r = Tj in the region r G 
[—T/2, T/2], and calculate the normalized determinant, using the following 
formula 0 for u±(T/2 ± Tj) 1. 

det'(-ff + V') </>(T/2) 

Aet(-dl+u%) XMT/2)’ { ’ ’ 

where the eigenvalues are determined with the use of the Dirichlet boundary 
conditions at ±T/2. The eigenvalue A is the one that appears in the valley 
equations and goes to zero in the limit T —> oo. The functions -0(r), ^ 0 ( r ) 
are the solutions of following differential equations: 

(-a? + V"(gMr) = 0, (B.2) 

{-d* + u= 0, (B.3) 


with boundary conditions, ijj(—T/2 ) = -0 O (—T/2) = 0 and ip(—T/2) = 
M~T/2) = 1. 

The function ijj is given by the following: 


V’oCO 



(^u+G+T/2) _ e -^+( r +' r / 2 )j ^ 


(B.4) 


which yields ^ 0 (T/2) ~ e UJ+T /(2 o;+). 

For the construction of the function we can utilize the auxiliary 

coordinate F(r), which is the eigenfunction with zero eigenvalue, 

(-d 2 T + V"(q)) F = 0. (B.5) 

with converging boundary conditions at r = ±oo. It has the following asymp¬ 
totic behavior: 

F(t) ~ F T e^ T ~ Tl \ (B.6) 
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for r —> ±T/2 as was shown in Section 2.2. We denote the other indepen¬ 
dent solution of Eq, (|B.5|) by F(t), in which we choose to have the following 
Wronskian: 

W = F(t)F(t) - F(t)F(t) = 2u+io-. (B.7) 

This implies that the asymptotic behavior of F(t) for r —> ±T/2 is the 
following: 

F(t) ~ F^e^ T - Tl) , (B.8) 

with F± = ^w^/F±. Using these solutions, the function is given as 
follows, 

< a9 > 

which satisfies the required boundary conditions at r = —T/2. 

Next we will evaluate the lowest eigenvalue A, 

(-d 2 T + U"(g)) $(t) = A$(t), (B.10 ) 

with ^(±T/2) = 0. We can construct the eigenfunction ^(t) as a perturba¬ 
tion series in A as follows, 


— /*+^V2 _ 

ip(r) = "0(t) T A / dr'Gir, 

J-T/2 

/*- 1 - ^ / 2 

= + A / dr'G(r,r')V’('T / ) + 0(A"), (B.ll) 

J-T/2 


where the function given by Eq.( [B.9| ) is used so that the boundary 
condition -0(—T/2) = 0 is satished. The Green function is defined 

by, 

G(t, t’) = A(f(t)F(/) - F(t')F(t))0(t - /). (B.12) 

and satishes the following, 

{-dl + V"(q)) G(t.t') = S(t - t‘). (B.13) 


The boundary condition -0(T/2) = 0 determines the eigenvalue A. Straight¬ 
forward calculation shows that the product of the second term of the Green 
function ( |B.12|) and the F{r') term in -0(r') makes the dominant contribution 
to ij(T/ 2). Therefore we find that 





a l/m r n A.-FV 

W2u+F(-T/2) J-r/2 { 


Combining the results obtained so far, we find that 
det \-dl + V") e G—^+){T/ 2 -n) -oo 

det(—+ w+) 2 U-F + F- J-o o 


(B.14) 


(B.15) 
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C The relation between the valley parameter 
and the separation of valley-instantons 

As has been shown by the explicit numerical integration of the valley equa¬ 
tions carried out in Section 2, the solution approaches to the well-separated 
pair of the (anti-)valley-instantons as the parameter A in the equation ap¬ 
proaches zero from below. In this appendix we will clarify the relation be¬ 
tween A and the separation R. For the sake of simplicity we will confine 
ourselves to the case of the II valley. Obtaining a similar result for the II 
valley is a simple matter. 

Assuming that A is small enough that the solution consists of a well- 
separated pair of valley-instanton and anti-valley-instanton, we can locate 
the valley-instanton at r = 0 and the anti-valley-instanton at r — R. At the 
symmetric point r 0 (= R/ 2), the time derivatives of the configuration vanish 
simultaneously; F(r 0 ) = q(r 0 ) = 0. In the asymptotic region to the far left of 
the valley-instanton, r < 0 and |r| 3> l/u;+, the solution is well described by 
linearizing the valley equations ( |2 .3 1|) - ([2 . 32|) . By requiring that it becomes 
a vacuum at q + as r —> — oo, the solution is described by two parameters F + 
and Q-, 

(C.l) 

(C.2) 



F + e K+ ' 

q-q+ - 

(q+- 


F± 

A 


F+ 

+ T -|_ —e K+T 

A 


where k± = \/oj± — A. Of two degrees of freedom in F + and Q + , one is related 
to the overall translation and the other should be fixed so that q and F are 
zero at To- We already know the asymptotic form of the valley-instanton as 
Eqs. (pAl)-(^60D. Thus 

Q+ F + -> -6 eg (C.3) 

9 

at the limit of A —> 0. 

Similarly the solution is well described by linearized equations around To, 
for the configuration is almost a vacuum at Since F(tq) = q(r 0 ) = 0, we 
can conclude the II valley behaves as 

F ~ F_e~ K ~ T + F_e K - ( ' T ~ 2T °' > 

( F - 
q-q- - [Q-- 


(C.4) 


A 


e~ UJ ~ T + —e~ K ~ T 

A 


Q_ _ ^Z.\ e u-{r-2r 0 ) + ^ eK _(r-2r 0 )_ 

A / A 


When A —> — 0, r 0 goes to the positive infinity and 


Q- -»• -- 

9 


F_ 


-6 eg, 


(C.5) 


(C.6) 
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obtained again by comparing (|C.4| ) and ( |C.5| ) with the asymptotic forms for 
positive r in (|2.59|) and ( |2.60|) . 

We use the conservation of Hy (|2.35|) to relate A to r 0 . At r = — oo Hy 
is equal to — V(q + ), while at r 0 

Hy ~ -V(q-) - l -u 2 _ [q(T 0 ) - q-f + jU 2 _F(t 0 ) [q(r 0 ) - g_] - ^F(t 0 ) 2 . 

(C.7) 

and (|C.5|) to relate 


Equating the two expressions for Hy and using 
F(t 0 ) and g(r 0 ) to F_, Q_, and r 0 , we obtain 


A = 


V(q+) - V(q _) L 


2u; 2 F_Q_ + F 2 (u-T 0 — 1) 


= -2 LU-TQ 


(C.8) 


Although F_ and Q- may have an 0(A) correction, they are exponentially 
small compared to the values in ( |C.6|) and are negligible. Thus we obtain 


A = 


24cj 2 + 36eg 2 (u;_-R - 2) 


.—id— R 


(C.9) 


The valley parameter A decreases simply as e 


—uj— R 


at e = 0, while it has an 


exotic dependence oj_Re UJ ~ R at non-zero e. 


D WKB evaluation 


In this appendix, we will re-derive the result ( |4.15|) by the WKB approxima¬ 
tion. The system we will consider obeys the Schrodinger equation 


1 d 2 

W (,) 


<h(g) = E®(q). 


(D.l) 


Around the upper minimum q + of the potential, this becomes 

Id 2 1 


2 dq 2 + 2 q 


<%) = E$(q), 


and a solution which vanishes at q —> — oo is 

<%) = AD U (-V2q), 


(D.2) 


(D.3) 


where E — v \ and A is a constant. The function D u is the parabolic 
cylinder function [^]. And around the lower minimum of the potential, 
the Schrodinger equation becomes 


H 2 1 

2 dq 2 




$(g) = E$(q), 


(D.4) 
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and a solution which vanishes at q —> oo is 

$(q)=AD v+e (V2(q~l/g)), (D.5) 

where A is a constant. We must now connect these solutions with that in 
the forbidden region. 

In the forbidden region, the usual semi-classical expression for the wave 


function is available: 

C\ r r q i C-2 r [i ] 

$(<?) = k( q y/2 exp [~ J q k ( x ) dx \ + k ( q y/ 2 exp [J q k K x ) dx \ > ( D -6) 

where 

k{q) = yj'q 2 (l - gq) 2 - 2egq - 2E . (D.7) 

and q t (i = 1,2) are the tunneling points V(qi) = E. Ci (i = 1,2) are 
constant. When eg 2 <C 1, these are given by 

gqi = eg 2 + \j2Eg 2 + (eg 2 ) 2 , (D.8) 

gq 2 = 1 + eg 2 - J2(E + e)g 2 + (eg 2 ) 2 . (D.9) 


The integral in the wave function can be evaluated by expanding about g 2 : 

rq 1 rgq . 

/ k(x)dx — — / [w 2 ( 1 — ta ) 2 — 2eg 2 ta — 2 g 2 E] lk2 dw 

dqi g 2 Jgqi 

1 /' 9 f / rz ew + E 

= — / tail — ta)dta — / — - --cha + • • • (D.10) 

g z Jgqi Jzi w( 1 — ra) 

1 rl 2 1 s / w \ o / 

= “2 o*" “ Q w - 9 Eln 1 - + re ln(l -w) + • • • 

g 2 L2 3 VI - wj \ gqi 

In the region near the upper minimum, it becomes 
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In the final expression, we have neglected the terms which vanish as g —> 0, 
so the wave function becomes 


$(<?) 


exp (-!,* + 

C 2? - 1/2 exp (±q 2 - E In -e]+- 

Cie-’ s / 2+E (2B)- E / 2 ? ‘' + C- 1 e‘ i ‘ l '‘- E {2E) E l' t q-‘' + ■ ■ •. 

(D.12) 


We are here using fc(g) = q + • • •. This wave function must be connected 
with the one given by Eq.( p.3| ). As the wave function ( p.3| ) has the following 
asymptotic form for q 1: 


<f>(g) = A 


=-9 2 /2 1 


-V2qY 




e ± ™ i e - q 2 / 2 (_ V 2 q )- u ~ 1 + 


T(- u) 

so we obtain the following matching condition 

±1^0-12-1/2 _ -(2i/+ 1)/ 9 , p-h+i/2) 

r (- i /) “ c 1 ! 1 + j 


, (D.13) 


(D.14) 


In a similar way, by connecting the wave function ( p.6| ) with the one 
given by (|D.5| ) near the lower minimum of the potential, we obtain another 
matching condition: 


_^_ e ±(v +e )iri 2 -(v+ 1/2+e) = 9l e ~^9 2 ( 2z/ _|_ 1) -(v+1/2) g -(4v+2+2e) . 

T(-u — e) C'2 


(D.15) 


Finally, if we eliminate the undetermined ration of the coefficients C\ /C' 2 in 
the matching conditions, the master equation ( |4 .1 5|) for the non-perturbative 
corrections of the energies is reproduced: 



(D.16) 


where 

p-1/63 2 

« =- TFT- 

gTT 


(D.17) 
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